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ABSTRACT 

Chandra high energy resolution observations have now been obtained from 
numerous non-peculiar O and early B stars. The observed X-ray emission line 
properties differ from pre-launch predictions, and the interpretations are still 
problematic. We present a straightforward analysis of a broad collection of OB 
stellar line profile data to search for morphological trends. X-ray line emission 
parameters and the spatial distributions of derived quantities are examined with 
respect to luminosity class. The X-ray source locations and their corresponding 
temperatures are extracted by using the He-like f /i line ratios and the H-like 
to He-like line ratios respectively. Our luminosity class study reveals line widths 
increasing with luminosity. Although the majority of the OB emission lines are 
found to be symmetric, with little central line displacement, there is evidence for 
small, but finite, blue-ward line-shifts that also increase with luminosity. The 
spatial X-ray temperature distributions indicate that the highest temperatures 
occur near the star and steadily decrease outward. This trend is most pronounced 
in the OB supergiants. For the lower density wind stars, both high and low X-ray 
source temperatures exist near the star. However, we find no evidence of any high 
temperature X-ray emission in the outer wind regions for any OB star. Since the 
temperature distributions are counter to basic shock model predictions, we call 
this the ^'near-star high-ion problem" for OB stars. By invoking the traditional 
OB stellar mass loss rates, we find a good correlation between the /ir-inferred 
radii and their associated X-ray continuum optical depth unity radii. We conclude 
by presenting some possible explanations to the X-ray source problems that have 
been revealed by this study. 
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Subject headings: stars: early-type - stars: X-rays - stars: winds, outflows - 
stars: shocks - stars: magnetic flelds - X-rays: stars 

1. Introduction 

The Chandra SateUite has provided astronomers the abihty to study the high energy 
resolution X-ray emission line spectra of numerous stars other than the Sun. With regards 
to early- type stars (hereafter OB stars), these line profiles allow us to study the complex 

stellar wind distribution of the X-ray source regions. The X-ray line profile shapes provide 
information on the Doppler velocities associated with the line formation regions, and allow 
us to comment on the shocks embedded in these stellar winds. In addition, emission line fiux 
ratios allow us to explore the radial distributions of the source regions and their associated 
X-ray temperatures. 

The main goal of this paper is to present the X-ray data from a large collection of 
"normal" OB stars and use the data to search for trends in the X-ray emission line parameters 
as a function of spectral type and luminosity class. Morphological trends in the X-ray 
emission line characteristics can help define important properties of X-ray emission that 
will be used to constrain the subsequent X-ray source modeling efforts. Although this type 
of study involves fewer assumptions than used in the modeling of each specific star, we 
suggest that it provides a broader, more comprehensive view of the X-ray source region 
characteristics. 

In our analysis we will also be interested in comparing how X-ray parameters that can 
be expressed in terms of velocities (i.e., line widths, line shifts, and X-ray temperatures) 
relate to the stellar wind terminal velocities and radial velocity structures. There are a 

variety of ways to discuss the data within the framework of several different shock pictures: 
1) an ambient stellar wind impinging on dense blobs (Lucy & White 1980); 2) a distribution 
of "saw-tooth" forward shocks (Lucy 1982); 3) a wind running into a driven wave which 
produces forward and reverse shocks (MacFarlane & Cassinelli 1989), and; 4) a detailed 
hydrodynamic numerical simulation of the line-driven instability which incorporates the 
complex time-dependence of the line driving force which also produces forward and reverse 
shocks (Owocki, Castor, & Rybicki 1988; Feldmeier 1995). The primary difference between 
these models is that the line-driven instability, or more appropriately the "de-shadowing 
instability", accelerates the pre-shock wind plasma to velocities that are larger than the 
ambient wind velocities resulting in a highly rarefied wind structure prior to the shock jump, 
and the post-shock velocity (i.e., the velocity of the X-ray emitting plasma) is predicted 
to be comparable to the ambient wind velocity. In the other shock models, the post-shock 
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velocity is some fraction (< 1, with a typical value of ~ 0.5) of the ambient wind velocity. 
Correspondingly, the post-shock density (i.e., the density of the X-ray emitting plasma) of 
the de-shadowing model is then also comparable to the ambient wind density, whereas, for 
the other shock models, the post-shock density is a factor of four times larger than the 
ambient wind density as determined by the Rankine-Hugoniot equations when the pre-shock 
Mach number is very large. Since our goal is to provide a global view of the OB stellar X-ray 
emission line properties and their dependencies on stellar parameters, we choose to use what 
we will refer to as a "basic shock model " description in our discussion of the relationships 
between the X-ray and wind velocity parameters. The premise of this model is in-line with 
the blob, saw-tooth, and driven wave shock models where the pre-shock velocity at any given 
radius is equal to the local ambient wind velocity as defined by the standard /3-law velocity 
structure. 

In §2 we provide a summary of the general OB stellar X-ray emission properties as deter- 
mined from previous studies. In §3 we present our selected stars, data reduction procedures, 
the X-ray lines that are used, and a description of the basic X-ray emission line parameters for 
our line fitting procedure. Using our Chandra High Energy Transmission Grating Spectrom- 
eter (HETGS) data along with the available archived Chandra data, a total of 17 OB stars 
are studied in this analysis. The observed velocity information (half-width-half-maximum 
and line-shifts of the line emission peak) displayed as histogram distributions, illustrating 
their dependence on luminosity class, are presented in §4. In §5 we present our discussion 
and analysis of X-ray emission line ratios and their application towards estimating the values 
of the source temperatures and the predominant stellar wind X-ray source locations. In §6 
we focus on deriving the radial temperature distributions of the OB X-ray sources. We find 
that this confirms the "near-star high-ion problem" found in earlier isolated stellar studies. 
The results and conclusions are summarized in 57. 



2. Summary of OB Stellar X-ray Emission 

Spectral lines and their profiles have been providing crucial information about OB stellar 
X-ray astronomy since the beginning of the field. Broad UV spectral lines of superionization 
stages such as O VI were discovered in Copernicus observations (Lamers & Morton 1976). 
These broad profiles and the persistence of the superionization stages over a well delimited 
range in spectral types led to the realization that these ions were produced by the Auger 
Effect in which two electrons are removed from the dominant stages of ionization following 
K-shell absorption of X-rays (CassineUi & Olson 1979). These authors suggested that the X- 
rays responsible for the superionization could come from a hot corona at the base of the cool 
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wind. Soon thereafter the Einstein Observatory discovered X-rays from OB stars (Harnden 
et al. 1979; Seward et al. 1979). However, the observations did not show the expected large 
attenuation of soft X-rays by the overlying cool wind. This led to the realization that the 
X-ray emission must arise from shock structures embedded within the stellar wind (Lucy & 
White 1980). The origin of these shocks is associated with the instability of line driven winds 
to velocity distribution disturbances (Lucy & Solomon 1970). In the Lucy & White case, the 
shocks are bow-shocks around radiative-driven clumps, whereas in the model of Lucy (1982), 
the X-rays are formed in a periodic shock structure. Spectral lines at X-ray wavelengths were 
detected by the Einstein Sold State Spectrometer (SSS) (CassineUi & Swank 1983). They 
found evidence of high energy ion line emission (Si xill & S xv) from the late 0-supergiant, 
C Ori. Since these lines are formed at such high temperatures, they suggested that these 
ions are located in magnetically confined regions at the base of the wind. In addition, 
since the overall X-ray fiuxes from OB stars are less variable than would be expected from 
spherical shells of shocked wind material, Cassinelli & Swank also suggested that the shocks 
were fragmentary in form, such that there would be a statistically steady supply of X-ray 
emission from these wind distributed "shock fragments". Some of these early ideas regarding 
X-ray source regions have persisted to the present time. 

The iiO^yir Satellite had a greater sensitivity than the Einstein Satellite, although not 
the energy resolution of the SSS instrument. In a survey of the Bright Star Catalogue OB 
stars, Berghofer et al. (1996 & 1997) confirmed that all 0-stars are X-rays sources, with 
X-ray luminosities following the "hot-star law" Lx/LboI — 10~^ until about Bl. For the 
later B stars, Cohen et al. (1996) found a rapid decrease in X-ray luminosity with spectral 
type which could be explained primarily by the slower speeds of the B star winds and the 
reduced mass loss rates. 

With Chandra, high spectral resolution of X-ray line emission from OB stars has now 
become observable with the HETGS Medium Energy Grating (MEG) and High Energy 
Grating (HEG) detectors. The impact of this diagnostic capability became very clear in the 
first detailed study of HETGS data from an 0-star ((^ Ori) by Waldron & Cassinelli (2001). 
Their analysis found three fundamental unexpected results associated with: 1) X-ray line 
profile shape characteristics; 2) the correlation between the wind X-ray source locations and 
their respective X-ray continuum optical depth unity radii, and; 3) the presence of deeply 
embedded high energy ionization stages. 

One of the biggest surprises from Chandra and XMM-Newton high energy resolution 
observations of OB stars is the absence of blue-shifted, asymmetric X-ray emission lines. 
These resolved spectral lines had been predicted by MacFarlane et al. (1991) to be broad 
and skewed towards the blue, owing to the fact that the long- ward ("red- ward") radiation 
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from the back side of the star is more strongly attenuated than the short- ward ("blue-ward") 
radiation from the near side. In the Waldron & Cassinelli (2001) analysis of ( Ori broad 
line profiles were seen, but these lines lacked the predicted skewness and blue-shifts. In the 
case of the 04f star ( Pup the lines arc somewhat similar to the predicted blue-shifts (Kahn 
et al. 2001; Cassinelli et al. 2001). However, nearly all other OB stars show minimal blue- 
shifted lines (within the resolving power of the instruments), and the lines are essentially 
symmetric (Waldron & Cassinelli 2002; Miller et al. 2002; Cohen et al. 2003; Mewe et al. 
2003; Schulz et al. 2003; Waldron et al. 2004; Gagne et al. 2005). Hence, explaining the 
lack of substantial blue-shifted lines that are presumably formed in an outflowing wind has 
become one of the most perplexing problems to emerge from high spectral resolution X-ray 
observations. 

Waldron & Cassinelli (2001) were the first to demonstrate that a stellar wind distribution 
of X-ray sources distributed above 1.5 stellar radii could explain the observed fines profile 
shapes seen in C Ori if the stellar wind mass loss rate is at least an order of magnitude 
smaller than previously thought. Since then, several studies (e.g., Kramer, Cohen & Owocki 
2003; Leutenegger et al. 2006; Cohen et al. 2006) have found that reduced mass loss rates 
are needed to explain the observed line profile shapes if the X-ray sources are originating 
from a distribution of stellar wind shocks (e.g., the shock model developed by Feldmeier 
1995), and the majority of the observed X-ray emission is found to arise from a stellar 
wind location between 1.5 to 2 stellar radn (Leutenegger et al. 2006; Cohen et al. 2006). 
At first this apparent mass loss rate reduction requirement was addressed as evidence that 
these winds are highly clumped, since a clumped wind would provide channels to allow the 
deeply embedded X-rays to escape more freely. Furthermore, a clumped wind would also be 
consistent with the mass loss determination results from IR and radio observations, since a 
lower mass loss rate could lead to the same emergent flux as an un-clumped wind of a higher 
mass loss rate. A key effect of wind clumping is to enhance the escape probability of X-ray 
photons, and this has also become known as the "porosity" of a stellar wind. However, as 
we will see later, for some of our stars, the reduced mass loss idea runs counter to a key 
X-ray observational result, the observed X-ray sources are located at their respective X-ray 
continuum optical depth unity wind radii as determined from using traditional mass loss 
rates (see discussions in Sec. 5.1 & 7). 

The perception that these stars do have reduced mass loss rates is based primarily 
on FUSE observations and analyses of the P V P-Cygni fine profiles from several 0-stars 
(FuUerton, Massa, & Prinja 2006). From this study, it now appears that either the mass loss 
rates are a factor of 10 or more lower than previously thought, or these winds are severely 
clumped over small spatial scales. Fullcrton et al. argue that UV lines from the ion P v 
(1118, 1128 A) should give line depths that are independent of clumping, consequently, their 



-6- 



results imply lower mass loss rates for these stars. However, if there are clumps in these 
winds and these are surrounded by X-ray producing shocks, it is questionable to assume 
that the ionization balance of phosphorus is not shifted to higher states by the Auger effect, 
or whether recombination in the high clump density decreases the ion abundance of P V. 
So the assumption that P V is dominant everywhere is arguable. With regards to X-rays, 
some authors find that clumping over small spatial scales appears to be able to explain 
the observed X-ray line profile shapes (Feldmeier, Oskinova, & Hamann 2003; Oskinova, 
Feldmeier, & Hamann 2004). However, Owocki & Cohen (2006) argue that the required 
porosity lengths are unhkely, hence, they conclude in favor of reduced mass loss rates. We 
find it ironic that 25 years ago the possibility of reduced mass loss rates would have permitted 
base coronae models to be an acceptable explanation of the X-ray emission from OB stars 
(Cassinelli et al. 1981; Waldron 1984). However, base coronae models were rejected because 
the needed mass loss rates disagreed with values derived from the available radio fluxes. 
Nevertheless, a better picture involving fragmented shocks in the winds emerged from the 
observations, and it is now abundantly clear that shocks are responsible for a major fraction 
of the X-ray emission from OB stars. However, there remain a large number of problems 
that we address in this paper that have led to questions about the overall nature of these 
wind distributed X-ray sources. 

The Chandra high energy resolution data has allowed us to exploit a major line emission 
diagnostic, the ratio of the forbidden to intercombination emission lines (//i) arising from 
He-hke ions (ranging from O Vll to S xv). This ratio had long been known as a very 
useful diagnostic tool for determining solar X-ray electron densities (Gabriel & Jordan 1969). 

However, it was clear that the densities derived with this interpretation turned out to be far 
too high (> 10^^ cm^^) for ^ Pup and ( Ori. Kahn et al. (2001) and Waldron & Cassinelli 
(2001) quickly realized that the cause of this weakening of the /—line relative to the i— line 
is not from collisional excitation, but rather radiative excitation from the presence of the 
strong UV/EUV flux from the photospheres of OB stars. The effect of radiation flelds on 
the f/i ratio had been accounted for in the calculations of Blumenthal, Drake, & Tucker 
(1972). The intense photospheric radiation near OB stars causes a depopulation of the 2^5"! 
level (/— hne) and a higher population of the 2^Pj levels (J = 0, 1, 2) (i— lines) by photo- 
excitation. This better understanding of the excitation process meant that the f/i ratio 
could be used to derive the radial distances of the predominant X-ray sources using the 
geometric dilution factor of the UV/EUV radiation. Basically, the smaller the f/i ratio, the 
closer the X-ray source is to the star. 

Waldron & Cassinelli (2001) found that the fir-inferred radii (Rfir) derived from analyses 
of the observed f/i ratios indicated that the He-like ions are present over a wide range 
in radial distances from the star. This means that there are X-ray sources distributed 
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throughout the stellar wind as would be expected from a distributed of stellar wind shocks. 
In addition, these Rfir were found to correlate with their respective X-ray continuum optical 
depth unity radii, Rr=i, i.e., the wind location where the X-ray continuum optical depth 
has a value approximately equal to unity. This observed correlation means that the highest 
energy He-like ions are located deep in the wind, with a steady progression outward of the 
lower energy He-like ions. This spatial distribution has now been observed in other stars with 
dense winds, C Pup (CassineUi et al. 2001), 5 Ori (Miller et al. 2002), and e Ori (Waldron 
2005). Somewhat surprisingly, this correlated behavior has also been observed in the highly 
luminous 0-star, Cyg 0B2 No. 8a, a star which is believed to have a mass loss rate at least 
5 times larger than ( Pup (Waldron et al. 2004). As discussed by Waldron et al. (2004), 
although the relationship between the Rfir and corresponding Rt=i is not exact, the general 
conformance is rather good. Our understanding of this behavior is related to the strength 
of the stellar wind opacity at a given energy. Since the opacity scales as ^ X^, at long 
wavelengths the wind can be very optically thick to X-rays. Hence, we can only see the long 
wavelength line radiation, i.e., the lower ion stages, that is emerging from the outer layers of 
the wind. In contrast, at short wavelengths, the wind becomes more transparent to X-rays 
and we can detect line emission that originates from deep within the wind. At these short 
wavelengths it is the higher ion stages that are producing the line emission. Since the X-ray 
cmissivity depends on the square of the electron density, even if the line emission for a given 
ion is arising from all depths in a wind, we would predominantly only see the emission from 
the highest density regions, which tend to be those as close to the star as possible. Because 
of attenuation by X-ray continuum opacity we expect to be able to detect radiation that is 
produced predominantly near optical depth unity. At this corresponding radial depth and 
farther out, the radiation can escape freely, but radiation from sources deeper in the wind is 
attenuated by the overlying material. The fact that there is agreement between Rfir derived 
from emission line ratios, and Rt=i derived from a consideration of the ambient wind optical 
depth, is certainly a relation that is to be expected. We ourselves had not predicted it, but 
found that it provided an explanation of the wide range of radii that are inferred from the 
observed f /i ratios. However, even this favorable relation has come into question since there 
is a call for the mass loss rates to be reduced by a large factor. We intend to see whether 
the Rfir and Rr=i arc correlated in this sample of many more stars than had been in our 
earlier studies of one or a few stars. 

This spatial distribution of He-like ions leads to what is perhaps the most significant 
problem in OB stellar X-ray astronomy. Very high He-like ionization stages such as Ar XVII 
(3.95 A) seen in Cyg 0B2 No. 8a (Waldron et al. 2004) and S XV (5.04 A) seen in C Pup 
(CassineUi et al. 2001) occur at wavelengths where the continuum opacity is very low, 
hence, any line emission from these ions can be detected from very deep within the wind. 
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The observed f /i ratios from these ions have confirmed that these high ion stages are in 
fact forming very close to the stellar surface (< 1.2 stellar radii). Prom the distribution of 
ionization stages, and the H-like to He-like line ratios to diagnose temperatures, Waldron 
(2005) found that the spatial distribution of X-ray temperatures (Tx) for OB supergiants 
steadily decreases outward from 20 MK near the surface to, 10 MK at 1.5 i?*, 5 MK at 3 
i?*, and 2.5 MK at around 8 i?*. Such a decreasing temperature distribution poses a problem 
in current X-ray studies because the temperatures required for the high ionization stages are 
higher than should be producible by shocks at such small radial heights in the wind! For shocks 
there is a maximal temperature which is determined by the jump in velocity across the shock 
front. One would expect that this velocity jump should be no more than the local wind speed, 
which is small near the base of the wind. Some authors (Leutenegger et al. 2006; Cohen et 
al. 2006) have presented arguments that this is not a problem because shocks can in fact 
form at the radii in question. However, that is not the full extent of the problem, not only 
are shocks needed, but the shocks must have the velocity jump sufficient to produce the hot 
temperatures and high ions that are observed to be originating near the star. Shock model 
predictions (e.g., Feldmeier, Puis, & Pauldrach 1997a; Runacres & Owocki 2002) show that 
shocks can form at and above 1.5 i?*, and the derived X-ray shock temperatures {Tx) at these 
low radial locations are found to be highly dependent on the line-driven instability triggering 
mechanism. For wind structures that are either self-excited or includes explicit photosphcric 
perturbations, the predicted shock Tx at these low radial locations are < 2 MK; well below 
the temperature required to produce the observed high ionization stages. However, as shown 
by Feldmeier et al. (1997a), photosphcric turbulence can generate a shock Tx of ~ 10 MK 
at > 1.5 i?*, but the contribution from this shock with regards to the overall observed X-ray 
emission appears to be quite small. Furthermore, although these models indicate that weak 
shocks can form below 1.5 i?*, these shocks cannot generate X-ray temperatures. 

Since the expected maximum shock temperatures are found to be significantly smaller 
than the temperatures required for these high ion stages, we call this the "near-star high- 
ion problem" (hereafter abbreviated as NSHIP). Possible explanations of the NSHIP 
have been proposed, but a consensus has not yet been reached. For the specific case of r Sco 
(BOV), Howk et al. (2000) suggested that clumps form from the density enhancements in the 
embedded wind shocks and these clumps become decoupled from the ambient wind velocity 
law. The clumps follow trajectories that have them fall back toward the star which could lead 
to high relative velocities between the clumps and the wind even at relatively small radial 
distances from the star. Each clump would also have a range of temperatures distributed over 
the frontal bow-shock extending from low Tx in the wings of the bow-shock to a maximum 
Tx at the frontal apex of the bow-shock. Howk et al. suggested that the conditions allowing 
for an in-fall of clump material might only be present in main sequence stars such as r Sco, 
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because wind clump drag forces would prevent the in-fall in more luminous stars with denser 
and faster winds. However, in light of the observational demand for clumpy winds, this 
model warrants a closer look. Very high Tx are seen in some OB stars such as 6^ Ori C 
(which has a field strength of several kilogauss; see Gagne et al. 2005) which are believed 
to have magnetically controlled outfiows analogous to Bp stars (ud-Doula & Owocki 2002). 
Even with more moderate fields there could also be solar-like phenomena occurring in these 
stars which may help explain the NSHIP. Studies have found that moderate magnetic 
field structures could rise buoyantly through the radiative envelope of OB stars (MacGregor 
& GassineUi 2003; MuUan & MacDonald 2005), and, in fact, the first complex magnetic 
topological map of the early B main sequence star, r Sco, has been revealed (Donati et al. 
2006). Hence, as had been proposed by GassineUi & Swank (1983), one could envision that 
these high ion stages may reside in magnetically confined loops close to the stellar surface. 
Alternatively there could be "coronal bullets" or "plasmoids" of fast moving material that 
might be ejected from the surface of the star owing perhaps to magnetic rcconncction in the 
sub-photospheric region as proposed for the sun by Cargill & Pneuman (1984). Although our 
study will primarily focus on the OB stars that do not show the extreme kilogauss magnetic 
fields with the hope of understanding the X-ray emission problems of ordinary OB stars, we 
include 9^ Ori G in our study so that comparisons can be made between "normal" stars and 
a highly magnetic one. 



3. Determining X-ray Emission Line Characteristics 

Existing detailed studies of individual OB stars have provided many interesting results 
about the stellar X-ray sources. However, a study of many stars, using identical analysis 
techniques is called for and is needed to search for general trends in the OB stellar X-ray 
source properties. In particular, we explore the line characteristics as a function of stellar 
luminosity class since stellar luminosity is considered to be the dominant contributor in the 
driving of these OB stellar winds. In this section, we present the OB stars used in our study, 
their relevant stellar parameters, and the X-ray emission lines used in our analysis. We also 
discuss our line fitting approach and the emission line parameters that can be extracted. 



3.1. The Stellair Sources and Observed Lines 

We have compiled the available archived HETGS data for 17 OB stars. Our analysis 
includes both the MEG and HEG spectral data. Although for most of our stars the HEG 
spectral lines have much lower signal-to-noise {S/N) than their respective MEG lines, we 
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use both data sets to check the consistency in derived parameters and possibly add support 
to the conclusions. Our program stars are listed in Table [1], along with the relevant adopted 
stellar parameters, and the Chandra observation identification numbers (Obs ID). For a few 
stars, their observations were carried out over two or more separate time segments. The 
MEG and HEG spectra and relevant spectral response files (ARFs & RMFs) were extracted 
using the standard CIAO software (version 3.2.2). For stars with multiple observations, these 
observations were co-added. Our study will focus only on the H-like and He-like lines and 
two Fe XVII lines. For the H-like and Fe XVII lines we consider only those lines with signal- 
to-noise ratio of S/N > 5. For the He-like fir (forbidden, intercombination, resonance) 
lines, we require that the total fiux from the three lines must have S/N > 5. In a few cases, 
primarily for the high energy lines, if a reasonable fiux has been established (i.e., > 3), we 
use these results only for estimating line ratios that provide interesting limits. The quantity 
listed in Tabled] is the scale factor associated with the stellar wind column density (see 
Table [1] notes for a definition). Table [2] lists the emission lines used, their rest wavelengths, 
and the temperature (T^) associated with the peak emission for each line. Also included 
in the table are the wavelength dependent X-ray continuum absorption cross sections for 
the stellar wind {<Jw), the ISM {(Tjsm), and several atomic parameters that are needed to 
calculate the He-like f /i ratio. The ajsM cross sections represent the "cold" gas limit, i.e., 
ISM absorption cross sections. However, it has long been known that for all O and early B 
stars, the value of aw is always less than ajsM for energies < 1.5 keV. This is because the 
intense UV/EUV radiation near the star increases the ionization state in the wind relative 
to the cold ISM. The difference in cross sections is especially large at low energies (e.g., see 
Fig. 2 in Waldron et al. 1998). While for energies above 1.5 keV (< 8.3 A), aw is ~ ctism- 
For example, from Table [2] we see that the difference between ajsM and aw first becomes 
less than 20% at wavelengths short-ward of 12 A(> 1 keV). The values of aw listed in Table 
|2] are representative of a typical 0-star at a T^jf = 35000/^. In general, aw at energies 
< 1.5 keV will be smaller for stars with larger Tg/j, and conversely for stars with lower 
Teff. For each star, the product Nwo x o'w gives the commonly used scaling parameter, r*. 
This optical depth parameter has often been used in various emission line profile modeling 
efforts (e.g., Owocki & Cohen 2001). To obtain the radial and wavelength dependent X-ray 
continuum optical depth, we consider a standard "/5-law" velocity structure with (3 = 0.8 
(Groenewegen, Lamers, & Pauldrach 1989). This gives r(r. A) = 5awWNwo{^ ~ w(r)°-^^), 
where w(r) is the radial wind velocity normalized by the terminal velocity (the procedure 
for determining Rt=i is discussed in the Appendix of Waldron et al. 2004). 



- 11 - 



3.2. Description of X-ray Emission Line Pcirameters 



We examine four basic X-ray emission line parameters: the total "observed "line flux, the 
"observed "line emission measure (EMx), the line width or more specifically the half-width- 
at-half-maximum [HW H M) , and for the line shift, we introduce a new terminology which 
we will refer to as the "peak hne shift velocity" (Vp) as discussed below. Line widths and hne 
shifts are expressed either in physical units (kms~^) or as normalized to the terminal velocity 
of the wind, v^o- In search studies for global trends, the important velocity parameters are 
those relative to the ambient wind, and not the actual physically velocities. We have to 
emphasize that the line fluxes and their associated EMx are the observed values since, in 
principle, the actual line fluxes and EMx of the X-ray emitting sources may be influenced 
by the presence of wind absorption. 

a) The total observed line flux is the total energy integrated line flux determined from 
our statistical best-fit modeling, and provides the line fluxes needed in our study of line ratio 
diagnostics. The extraction of a line flux for a given wavelength region must be handled with 
care if one wishes to compare various observed line flux ratios with theoretical predictions. In 
particular, we find that it is extremely important to account for all of the lines that could be 
contributing to a given wavelength sector of the line, such as other strong lines and satellite 
hues. The importance of including other hnes was clearly demonstrated by Ness et al. (2003) 
in their analysis of the He-hke Ne ix fir lines from Capella. Also, appropriate adjustments 
for contaminating lines must be made in the model line flux ratios before comparisons with 
observed ratios can be made as discussed in Section 5. 

b) In general astrophysical studies, if distance and the temperature (which provides a 
line emissivity) are known, the total line flux (Fx) can be used to extract an "observed" line 
emission measure using 



where d is the stellar distance, and e(T£,) is the maximum line emissivity (e.g., see Kahn 
et al. 2001) evaluated at Tl. However, in the case of OB stars, the interpretation of the 
total observed line flux can be somewhat ambiguous, especially for OB stars with massive 
winds, because there can be unobservable sources located deep within the wind. To account 
for the total amount of hot material in a wind, one would need to find the attenuation 
of the line flux by the overlying wind. In addition, this definition of EMx assumes that 
all the emission for a given line is at the same temperature throughout the wind, and as 
discussed in Section 6, we now know that there is a radial dependent X-ray temperature 
structure. However, these two effects would require a significant amount of modeling of the 
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wind structure, and the introduction of significant uncertainties. In this paper we want to 
focus on a uniform discussion of observational properties. Therefore, we choose to deal only 
with the "observable" EMx- Of course one who is interested in the total X-ray productivity 
of the wind needs to understand that our derived EMx are to be strictly treated as lower 
limits. 

c) The breadth of a line is expressed as a HWHM. This provides information regarding 
the dispersion in the velocity of the line emitting plasma. For the general case in which there 
is a spatial distribution of X-ray sources, the derived HWHM represent integrals over both 
depth and impact parameters within each small range in the linc-of-sight (LOS) velocity. 
However, even in the absence of details regarding the line formation process, we find useful 
information can be derived from the line widths. For example, if a line is only produced 
very close to the star, one would expect it to be narrow, with a width perhaps comparable 
to the expected thermal or turbulent speed at the base of the wind. In fact, this is often 
seen in the X-ray observations of cool stars such as Capella. In an OB star, a narrow line 
could also result from certain asymmetric wind structures such as X-ray material confined 
to volume sectors which are seen inclined relative to the observer's LOS. Even if the total 
outfiow velocity in these sectors may be large, the line width could be highly dependent 
on the inclination angle. For example, an observer's LOS perpendicular to the flow would 
essentially see a near zero line width. We find that almost all OB stellar X-ray emission 
fines have moderate to large HWHM, although a few stars do have very narrow lines (e.g., 
T Sco; Cohen et al. 2003). 

d) We have introduced the phrase, "peak-line shift velocity" (Vp) to represent the re- 
quired Gaussian line profile model velocity-shift needed to obtain a best-fit to the given 
observed line profile. We use this terminology to emphasize that, in general, Vp does not 
correspond to a Doppler shift of any specific part of the wind or atmosphere. The observed 
shift is affected both by the spatial distribution of the X-ray sources in the wind, and by 
the degree to which the wind absorption attenuates the radiation from each X-ray source 
region. The observed peak-line shift is weighted by the dominant source locations, from 
which the emitted X-ray line radiation can escape through the overlying wind. Consider the 
simplest case of a single spherically symmetric shell of X-ray emitting material moving at 
some velocity V . For an optically thin wind, the observed line profile would be symmetric 
and fiat topped extending from about —V to -\-V (neglecting stellar occultation effects), 
with negligible fine shift, Vp ~ 0, and a HWHM < V. However, for the case of an optically 
thick wind, the red- ward side of the line is more heavily attenuated, and thus the peak in the 
emission, Vp, would occur at a blue- ward shift. This velocity could be slightly smaller than 
V , and the line shape would be sloped down long- ward from this Vp, producing a triangular 
shaped line as had been predicted by MacFarlane et al. (1991). For this thick wind case, the 
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HWHM can be significantly less than V. For even more complicated scenarios, such as a 
stellar wind with many discrete X-ray source regions, the observed line shift Vp represents an 
"average" over all the source regions that are both capable of contributing to the observed 
line emission which are located at an X-ray continuum optical depth of about unity or less. 

3.3. Line Fitting Procedure 

To obtain a totally unbiased collection of emission line parameters we want to use a 
model-independent extraction method. Two methods have been considered: 1) Gaussian 
fits to the line profiles, and; 2) the "moment" method discussed by Cohen et al. (2006). 
In the latter, one calculates the first three moments of a line profile using only the actual 
observed count spectrum (no ARF and RMF corrections) which in turn provide information 
on the line shift, line width, and asymmetry of the line profile. Although this would appear 
to be a good unbiased way to describe a line, there are limitations as discussed by Cohen 
et al. (2006): 1) with regards to obtaining the HWHM from the 2"*^ moment, the moment 
analysis does not allow one to separate the effects of instrumental broadening from physical 
broadening, and; 2) the moment method cannot extract reliable information from blended 
lines (e.g., the He-like fir lines). Furthermore, the moment method cannot be used to 
extract the total line fiux. Since our goal is to provide a collection of physical line emission 
parameters and line fiuxes for single and blended lines in a self-consistent manner, wc chose 
to use the Gaussian line fitting procedure described by Waldron et al. (2004). However, we 
do find that the 1** moment method, which is used to extract Vp, provides almost identical 
results as compared to those obtained using the Gaussian method. 

For the Gaussian line fitting procedure, wc assume that all lines within a given wave- 
length region have Gaussian line profiles superimposed on a bremsstrahlung continuum. We 
use statistics to determine the best fit parameters (line flux, HWHM, Vp, & EMx)- We 
also considered the C-Stat (Cash 1979) approach to line fitting but did not find any advan- 
tage over the approach used in this analysis, basically because the C-Stat method was 
developed for handling weak lines and in this paper we are considering only the strongest 
lines. The error bars for each parameter are determined from the 90% confidence regions. 
For all line fits we assume a continuum temperature of 10 MK. The actual value assumed 
is not critical to our results since we are only interested in fitting the shape of the line and 
the strength of the line emission relative to the continuum. All line flux emissivities and 
line rest wavelengths are taken from the Astrophysical Plasma Emission Database (APED; 
Smith & Brickhouse 2000; Smith et al. 2001). It is worth stressing again that even though 
we have only chosen to analyze the strongest observed lines, many of the line regions contain 
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overlaps with other fairly strong hnes. We find it important to include all extraneous lines 
that may be contributing in a given wavelength region in order to get an accurate fit to the 
line profile of interest. In the case of the three He-like fir lines, all three lines are likely 
to be formed under the same physical conditions, and our model fitting procedure assumes 
that the HWHM and V-p are the same for all three lines, and only the line strengths are 
different. The model count spectrum is determined for the default wavelength binning of the 
MEG (0.005 A) and HEG (0.0025 A), using the appropriate ±1"* order ARFs and RMFs. 
To obtain the best fit model parameters and associated errors, both the model and MEG 
and HEG ±1** order spectra are re-binned to bin sizes of 0.02 A and 0.01 A, which are 
respectively the approximate resolution limit of the MEG (0.023 A) and HEG (0.012 A). 
The main advantage of this approach is that it provides tighter constraints on the model 
parameters. The determination of line ratios (discussed in §6) are based on the resultant 
best-fit model line fiuxes. 



4. Distributions Versus Luminosity for the HWHM, Vp, and EMy^ 

We provide histograms of the HWHM and Vp in physical units (kms~^) for the case of 
comparing all OB stars, and HWHM and Vp normalized to their respective star's terminal 
velocity (voo) in our study of luminosity class dependence. It is these normalized distributions 
that are of interest for developing an understanding of the behavior of the X-ray emission line 
parameters versus OB spectral types and luminosity classes, but is not necessarily relevant 
when looking at the overall OB distributions for all classes. The EMx histograms represent 
the physical observed emission measure (cm^^). 

First we show the line parameter distributions for all 17 OB stars, then we examine how 
both the MEG and HEG derived parameter distributions depend on the stellar luminosity 
class. We consider three luminosity class groupings: 1) main sequence, MS (luminosity class 
V); 2) giants (luminosity classes IV & HI), and; 3) supergiants (luminosity classes II & I). For 
all histograms, the bins represent the percentage of lines within a given range of HWHM, 
Vp, or EMx- The HWHM and Vp bin sizes are 200 kms~^ for the actual velocity values, 
and 0.1 for the normalized case. For the Vp histograms, the bin spacing is set up to center 
on Vp = 0. 

The observed EMx are derived from the total line fiux by using the approach discussed 
by Kahn et al. (2001), using the APED emissivities and eq. [H The major assumption used 
in extracting EMx is that we have assumed each line is at its maximum emissivity level 
(i.e., maximum temperature Tl). In our study of X-ray temperatures (see Sec. 5.2) we find 
that this is not a bad assumption since all extracted temperatures are very close to their 
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maximum values. We present the EMx histogram distributions in units of the log EMx in 
increments of 0.5. 



4.1. HWHM, Vp, and EMx Distributions for All OB Stars 

The HWHM, Vp, and EMx MEG and HEG histograms for all OB stars are shown in 
Figure [H Several key features are noted: 1) the MEG and HEG Vp distributions are nearly 
symmetric (both show slight asymmetry blue-ward) around Vp = with ~ 80% of all lines 
lying between ±250 kms~^; 2) the MEG and HEG EMx distributions indicate that more 
than half of all lines lie within a relatively small range of ~ 1.5 dex in log EMx] 3) the 
MEG and HEG HWHM distributions show a large range from to 1800 kms~^, and; 4) 
the MEG HWHM shows a double peaked distribution with peaks at 350 kms~^ and 850 
kms~^ which is not seen in the HEG distribution which only shows a peak at 350 kms^^. 
The explanation is related to the sensitivity of the instruments in that the MEG 850 kms~^ 
peak is due primarily to lines from low ion stages which are inaccessible or very weak in 
HEG spectra. 

The Vp distributions illustrate one of the most surprising results arising from HETGS 
data analyses, the majority of all OB lines show essentially no line shifts. We emphasize 
that almost all Vp are typically within the wavelength resolution limits of the MEG (0.024 
A) and HEG (0.012 A). For example, the MEG velocity resolutions at wavelengths of 25, 20, 
15, 10, and 5 A are respectively 276, 345, 460, 690, and 1380 kms^^ Correspondingly, the 
HEG limits can be obtained by taking half of the MEG determined values (long-ward of 20 
A is out of the HEG wavelength range). 

These OB stellar distributions illustrate a basic dilemma in early attempts to understand 
the X-ray emission from OB stars. As in pre-launch expectations (e.g., MacFarlane et al. 
1991), the expected hne broadness is observed, but there is clearly a lack of substantial blue- 
shifted line profiles. From this sample, it is now clear that this non-shifted line behavior 
holds for essentially all OB stars. Hence there must be some common property related 
to X-ray production that needs to be determined. In the following subsections we explore 
the HWHM, Vp, and EMx distributions versus luminosity class to search for a better 
understanding of the non-shifted behavior seen in the distributions of Figure [H 
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4.2. HWHM Dependence on Luminosity Class 

The MEG and HEG histograms showing the HWHM dependence on the three lumi- 
nosity class groups are shown in Figure [2] (normalized to Voo)- With regards to all luminosity 
classes, a key finding is that all HWHM are < Voo- We see that the histogram distribu- 
tions are similar for supergiant and giants with peaks between 0.35 to 0.45, and both show 
asymmetries towards lower values (except in their HEG distributions). The most notable 
difference is the HWHM /vod histogram for the MS stars where there is no sharp peak. In- 
stead, a rather a large range of 0.1 to 0.5 x Voo is evident. A partial explanation of this can 
be seen from an inspection of the observed f oo and N]yo listed in Table [TJ We see that the 
supergiant stars all have similar Voo and N]yo, but the MS stars have a rather large range 
in Voo- For the MS stars the broad peak may be due to the fact that their Nwo are weak, 
and thus the wind absorption effects do not play a major role in determining the line profile 
shapes. The giants also show a larger range in both Voo and Nwo- 

The actual velocities from the HWHM values are consistent with X-rays arising from 
shocks located in the accelerating parts of the wind and/or in the case of rotationally dis- 
torted fiows, from sectors of the wind with low velocities along the line of sight. Surprisingly, 
both the supergiants and giants show significant line emission at HWHM /voo < 0.3. At 
least for the supergiants, we would have expected to see very little emission at these values 
since, due to their dense winds, one simply cannot see to the base of the wind owing to the 
strong wind attenuation. A possible explanation is that if one cannot see to the back side 
of the star owing to wind attenuation, this would also tend to make the HWHM smaller, 
but this would also require a large blue-shift in Vp, which is not observed. For the MS stars, 
although they also have a large range in v^o, their HWHM /voo histogram is not affected. 
It may be that the low mass loss rates of the MS stars relative to the giant and supergaint 
stars leads to lower wind column densities, hence, the observer can see deeper into the MS 
winds. Thus we are primarily detecting X-rays from regions of low velocities instead of re- 
gions with speeds near Voo- Whereas for the supergaint and giant stars, depending on the 
line wavelength and the X-ray opacity, we can observe to a broader range of depths and have 
a HWHM that depends on Voo- Furthermore, since none of the groups show significant 
number of lines with HWHM/voo > 0.5, the majority of the observed X-ray emission must 
arise relatively close to the star (< 2 stellar radii) which means that strong wind shocks in 
the outer wind regions do not exist as predicted by early shock models. 
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4.3. Peak Line-Shift Dependence on Luminosity Class 

The MEG and HEG histograms showing the peak line-shift, Vp, dependence on the 
three luminosity class groups are shown in Figure [3] (normalized to foo)- Most notable, 
with respect to the maximum of the Vp distribution, is that each luminosity group has 
a maximum that is nearly identical to the sample as a whole. That is, each class has 
a maximum in its distribution occurring at very small velocities, which are less than the 
MEG and HEG spectral resolutions. Next, in regards to the skewness or asymmetry of the 
distributions, there is an increasing blue-ward asymmetry of Vp with luminosity group. This 
indicates that a small fraction of the lines do indeed exhibit finite blue-shifts. The blue-ward 
asymmetry extends only to about -0.15 Voo for MS stars and -0.20 Voo for the giants. For 
the supergiants the Vp asymmetry distribution is well pronounced up to approximately -0.35 
foo- Nevertheless, this is still far less than had been expected from early line calculation 
modeling by MacFarlane et al. (1991), and from the empirical shock model calculations of 
Owocki & Cohen (2001). 

The observed peak-shifts can be deduced from the uniform wind modeling expectations 
in several ways: 

• If the wind mass loss rate is less than the value inferred from radio observations (as 
suggested by Fullerton et al. 2006), this reduction in wind density would allow more 
X-ray radiation to emerge from the back side of the star (i.e., the red-shifted emission). 
This approach was first demonstrated by Waldron & Cassinelli (2001) in their analysis 
of C Ori, and used by Kramer et al. (2003) to model the X-ray line profiles and peak- 
shifts from ( Pup. The Cohen et al. (2006) re-analysis of ( Ori has also confirmed 
that this approach can provide a fit to these lines. However, as mentioned before, there 
are some questions about the new M estimates associated with the actual ionization 
fractions of phosphorus. At the present time we choose to adopt the traditional mass 
loss rates until the major change in mass loss properties has been verified and widely 
accepted. 

• If the X-rays can escape more easily from the sides of a shocked region, then the X-ray 
plasma produced along a LOS perpendicular to the observer (i.e., from the wind on both 
sides of the star) would dominate the emergent X-ray emission. This is the solution 
offered by Ignace & Gayley (2002), who used Sobolev escape probability theory to 
derive the line profile shapes. However, it is questionable that the Sobolev theory they 
used is valid for non-monotonic cases with discontinuities in the velocity, and jumps 
in the velocity gradient. Furthermore, the X-ray formation region is almost surely not 
being accelerated at the same rate that the stellar UV radiation is accelerating the 
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wind that is colliding with the shocks. Perhaps the treatment could be improved using 
the Rybicki & Hummer (1978) approach for non-monotonic velocities. Nevertheless, 
the results of Ignace & Gayley do illustrate that by having the line radiation escape 
out from the sides of a shock provides a plausible explanation of the symmetric line 
problem. 

• If the winds are clumpy, photon mean free paths can be increased and one can see more 
easily to the back side of the star. This increased porosity effect has been investigated 
by Feldmeier, Shlosman, & Hamann (2002). It also addresses the reason why the radio 
fluxes could be larger than inferred from a laminar wind case. Although, it appears 
that one should expect to sec more fiat topped X-ray emission lines than arc actually 
observed from this initial picture, Oskinova, Feldmeier, & Hamann (2006) have made 
improvements which appear to explain the minimal blue-shifts, provided that these 
winds are highly clumped. 

• If the outflowing stellar winds are geometrically conflned, then the observed line profllcs 
will be dependent on the observer's orientation relative to the geometric structure of 
the wind. This approach was studied by MuUan & Waldron (2006) by considering a 
two-component wind structure (i.e., a polar and an equatorial wind), where the polar 
wind is slower and less dense than the equatorial wind owing to the fact that the 
polar wind is hindered by surface magnetic structures. Their results indicated that the 
observed line-shifts will be dependent on the LOS where a pole-on view would yield 
minimal blue-shifts. The main advantage of this model is that a fairly large sector of 
the outflowing wind does not have to be clumped. 

There are surely other more complicated scenarios that could explain the properties of 
the observed line shifts, but these suggestions cover the recently studied ideas. Although 
several of these can explain the observed lack of substantial blue-ward peak shifts, a con- 
sensus has not been reached. Furthermore, this issue of un-shifted hues became even more 
problematic with the HETGS X-ray emission hne analysis for Cyg 0B2 No. 8a (Waldron et 
al. 2004). This star is beheved to have an M that is at least 5 times larger than any other 
previously studied OB star, which implied that large blue-ward peak shifts should be seen. 
However, the observations show that the X-ray lines from Cyg 0B2 No. 8a are similar to 
other OB stars in that only minimal peak shifts were observed. A particularly perplexing 
question is why do substantial blue-ward peak shifts show up only in Pup? This star has 
for many years been the prototypical early 0-star in regards to its wind properties, and 
observations of it played a major role in the development of line driven wind theory. 
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4.4. EMx Dependence on Luminosity Class 

The MEG and HEG histograms showing the EMx dependence on the three luminosity 
class groups are shown in Figure HI There are several interesting results shown in these 
distributions. First we should recall that EMx is a measure of the X-ray density squared 
times a volume element, and EMx scales as (M /vooY /R*- With this in mind, the results 
shown in Figure H] suggest the following, using the data listed in Table [H 1) Since the 
supergiants have very similar M, Voo-, and -R* then they should have almost identical EMx 
and indeed that is exactly what we see (the exceptions are Cyg 0B2 Nos. 8A & 9). 2) The 
giants have the largest ranges in M and Voo hence, what we see in their EMx distribution 
is clearly consistent with this spread in wind parameters. 3) The MS stars show a double 
peaked distribution. The explanation is clear by inspection of their wind parameters, which 
shows two extreme groupings of M and Voo (i-e., early MS stars have larger M as compared 
to late MS stars). In principle, these EMx values which were determined solely by fitting 
the observed line profiles, independent of any wind parameters, should provide a means for 
establishing X-ray source densities if we know the radial location, temperature, and geometric 
extent of the region (e.g., cooling length). This is beyond the scope of this paper but we 
plan to explore this possibility in a later paper. 

5. Line Ratio Diagnostics 

In astrophysical studies line emission ratios are commonly used as diagnostic tools. In 
our study of high spectral resolution X-ray astronomy, we focus on two line emission ratios 
that will be used to establish X-ray source spatial locations and temperatures. Spatial 
information is obtained from the ratio of the He-like ion forbidden to intercombination 
{f /i) emission lines. This ratio provides a diagnostic for estimating either the distance 
of the dominant X-ray emission zone from the central EUV/UV radiation source (radiation 
dominated case) or the electron density of the X-ray emission region (coUisional dominated 
case). Although we commonly refer to the z— line as if it were a single line, it actually consists 
of two lines that are unresolvable in both the MEG and HEG. All theoretical f /i line ratio 
calculations include the total emission from both lines. To determine X-ray temperatures, 
we use the temperature sensitive line ratios from H-like ions to He-like ions (here abbreviated 
as the H/ He ratio). The H/He line ratio increases dramatically with temperature as shown 
in Figure [51 As discussed by Waldron et al. (2004), the temperatures derived from H/He 
ratios provide a measure of the "average" temperature of the He-like ions and the more 
highly ionized H-like ions. 

Another temperature sensitive line ratio is the He-like G-ratio defined as [i + /)/t. 
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This ratio (introduced by Gabriel & Jordan 1969) has been used extensively as a temperature 
diagnostic in solar X-ray studies, and has been used in early studies of 0-stars (e.g., Schulz 
et al. 2000; Waldron & Cassinelli 2001). The G-ratio dependence on temperature is opposite 
to that of the H/ He ratio, the G-ratio decreases with increasing temperature and the slope 
versus temperature is much weaker than is the case for the H/He ratio. An advantage 
of the G-ratio is that it provides a temperature of only one ion, the He-like ion, but, as 
shown by Waldron et al. (2004), significant differences in the derived G-ratio temperatures 
as compared to the H/ He derived temperatures and their associated were found. These 
discrepancies may be related to "line-blending" effects and/or resonance line scattering in 
the r-line as discussed by Porquet et al. (2001). Complications in interpreting G-ratio 
derived temperatures have also been discussed in studies of late-type stars (e.g., Ness et 
al. 2003). Recently, Leutenegger et al. (2007) state that they have found evidence for 
resonance line scattering among the low energy He-like ions in their analysis of the C, Pup 
XMM — Newton RGS spectra. Although we provide a tabulation of the observed MEG 
and HEG G-ratios (see Table [7]), we believe that until we have a clear understanding of the 
G-ratio idiosyncrasies, it is premature at this time to tabulate derived G-ratio temperatures 
as it may lead to confusion and misinterpretations. 

Line ratio diagnostics require good energy resolution of the observed lines. Even with 
the high energy resolution capabilities of the HETGS, we still must allow for possible con- 
tamination from the blended dielectronic satellite lines and any other lines that are within 
the instrumental resolution limits (these lines are easily identified in the APED data tables). 
These line-blending effects are most pronounced when dealing with ratios that use at least 
one or more of the He-like fir lines. For H/He line ratios, the effects of line-blending are 
minimal, with the largest effect occurring for neon. The line-blending effects on the f /i line 
ratio are also expected, but will not be explored in this paper. The main reason is that 
the f /i ratio is dependent on three parameters (UV/EUV flux, spatial location, and X-ray 
temperature), whereas the H/He ratios are essentially only dependent on one parameter, 
the X-ray temperature. This allows us to easily tabulate the expected line ratios, including 
line-blending, that we use when comparing with observed H/ He ratios. Although the H/He 
ratios may also be somewhat dependent on wind absorption and resonance line scattering 
effects (see Appendix), all H/He ratio derived temperatures {TuHe) presented in this paper 
only include line-blending effects. 
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5.1. Analysis of He-Like f ji Ratios 

The most widely used line ratio diagnostic to have emerged from Chandra studies of OB 
stellar X-ray emission is the He- like f /i line ratio. This ratio has proven to be a valuable 
observable because, for the first time, we have a direct means of finding the predominant 
X-ray source stellar wind locations. As we have seen. X-rays can arise from a wide range 
of wind radii. However, as is the case in line formation in stellar atmospheres, there is 
a depth for which the contribution function is maximal. Since the X-ray line emission 
varies as the square of the density, deeper layers tend to contribute more strongly to the 
observable line strengths. It is information about the predominant zone in radial distance 
that we can obtain information from using the f /i ratio. Ever since the beginning of OB 
stellar X-ray astronomy, determining the location of this X-ray emitting plasma has been 
considered to be the most crucial quantity required for understanding the X-ray emission 
process in OB stars. The Einstein and RO S AT saXeWiies provided only low resolution data, 
hence, the only information regarding the source location was estimated from the change 
in X-ray attenuation across the major continuum jumps, such as the Oxygen K-shell edge. 
Historically, this was sufficient to show that the X-rays from OB supergiants could not all 
be coming from a thin coronal zone at the base of the winds (Cassinelli et al. 1981). 

The He-like f /i line ratios have been used in solar coronal studies as a diagnostic of the 
electron density of the X-ray emitting medium as first demonstrated by Gabriel & Jordan 
(1969). As the electron density is increases, a larger fraction of the ion is excited from the 
2^S level (the ground state of the triplet sequence) to the 2^Pj levels which results in a 
diminution of the strength of the 2^5* l^S transition (/—line) and an increase in the 
strength of the 2^Pj l^S transition (i— line). As discussed previously, the total z— line 
emission actually consists of two transition from the 2^Pj levels (J = 1 & 2). Although the 
J = level does not contribute to the i— line emission since it is strictly forbidden, it is a 
contributing transition in the rate equations which establish the relative populations of these 
levels (e.g., Gabriel & Jordan 1969; Blumenthal et al. 1972). 

In all O stars and early B stars the enhancement of the line relative to the /—line is 
caused by photo-excitation from the 2^Si level by the strong UV/EUV photospheric emis- 
sion. The radiation required depends on the ion, and ranges from 1638 A (O vil) to 674 
A (S xv) (see Table [2]). The photospheric lines required to excite O vil Ne ix and Mg xi 
all have wavelengths long- ward of the Lyman limit (directly observable regions), whereas 
the lines required to excite Si xill and S xv have their wavelengths in the EUV, short-ward 
of the Lyman limit (unobservable). Hence, theoretical f /i ratios for Si xill and S xv are 
dependent on model atmosphere fiuxes. If one of these excitation line wavelengths happens 
to coincide with a photospheric absorption line that for some reason is deeper than predicted 
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by the photospheric models, the distance inferred from the f /i ratio would be even closer. 
Conversely if the overall EUV continuum is larger than expected from model atmospheres the 
formation region is located farther out. An example of these effects is shown by Leutenegger 
et al. (2006). 

In the early Chandra studies by Waldron & Cassinelli (2001), Cassinelli et al. (2001), 
and Miller et al. (2002), we felt assured that the we were not too far off in our radial distance 
estimate by the fact that the radii derived from f /i (Rfir) and from the X-ray continuum 
optical depth unity radii (i?r=i) turned out to be nearly the same. However, it is important 
to point out that this correspondence is not exactly one-to-one, as discussed by Waldron et 
al.(2004). Furthermore, the recent arguments by FuUerton et al. (2006), and others, that 
the mass loss rates of OB stars are incorrect, are a source of concern and confusion as to 
why the Rfir and Rr=i are in reasonably good agreement. We will return to this later. 

Following the approach used by Blumenthal et al. (1972), the radial dependence of the 
f/i ratio (commonly labeled as R) can be determined by 

R(r] = — l'2'l 

^ ' 1 + n,ir)/Nc + W{t)^I<^c ^ ^ 

where is the photo-excitation rate from 2^5*1 2?Pj\ 

(3) 



and Rq represents the low density and = limit, ne(r) is the X-ray electron density of the 
X-ray emitting plasma, Nq is the critical density, (pc is the critical photo-excitation rate, 
and Ui, is defined as the surface photospheric radiation flux density which allows us to factor 
out the radial dependent geometric dilution factor, W{r) = 0.5(1 — (1 — (R^/rYY^'^). The 
values of Nc, (pc, and wavelengths of the 2^ Si — > 2^Sj transitions are given in Table [2l 
Although 0(7 is only dependent on atomic parameters, Nc is inversely proportional to the 
collision rate. Hence, Nc is also weakly dependent on the temperature, and the values listed 
in Table [2] are evaluated at each ion's expected maximum line emission X-ray temperature 
(Tl). 

The format of eq. [2] allows us to see the distinction between coUisional domination and 
radiation domination which is determined by the relative strengths of 0/0c and ne{r)/Nc- 
For OB stars, the radiation term is dominant throughout the wind except when one is 
interested in conditions very close to the photosphere. Assuming that the radiation term 
is dominant, we can neglect the density term and obtain a relationship between radius, the 
stellar photospheric radiation (i.e., the stellar effective temperature, T^jf), and the observed 
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f /i ratio (R), which is given by 



Wir) . ^(^ - 1) ,4) 



The main advantage of this equation is that for any OB star with an observed f /i ratio we 
can estimate a radial distance corresponding to the source of He-hke ion emission, provided 
that we know the corresponding radiative fluxes at the three A/_j (see Table E]) which can 
be determined using stellar atmospheric model spectra (e.g., Kurucz 1993; Hubeny & Lanz 
1995). Leutenegger et al. (2006) found that there are model atmosphere dependent differ- 
ences in the f /i ratio. We find that these differences appear to be relatively minor and we 
chose to use the Kurucz (1993) model atmospheres in our calculations. 

The f /i ratio defined in eq. [2] is referred to here as the ^^localized^^ f /i ratio. To 
understand the difference between this localized f /i ratio approach and the more complicated 
scenario, consider a single small isolated test packet of X-ray emitting plasma, at some 
constant X-ray temperature (Tx), which can be placed at any wind radial location. As 
this isolated test packet moves outward through the wind towards larger radii, the f /i ratio 
increases because there is now less depletion of the /—line and less enhancement of the 
i— line. This change is due to the radial decreasing strength of the radiation field (i.e., the 
decreasing dilution factor). In general, the f /i ratio is also weakly dependent on Tx through 
the critical density term since this term is inversely proportional to the collision rate, and 
obviously, the strengths of the fir lines are dependent on the fractional abundance of the 
He-like ion as determined by the temperature of the X-ray emitting plasma. 

Now consider a more general f /i case, the distributed^^ f /i where there is a radial 
distribution of X-ray sources throughout the wind. Here, as compared to the localized case, 
the total line fluxes of both the /— and z— lines are determined by an integration process (see 
Leutenegger et al. 2006), or equivalently, a summation of a very large (or infinite) number of 
isolated test packets distributed throughout the wind. However, since the strength of each 
line emissivity scales as the packet's density-squared, the contributions to /—line deep in the 
wind are depleted by radiative excitation, and only the i— line accumulation benefits from 
the density-squared effect since the i— line is much stronger than the /—line at low radial 
distances. On the other hand, the flux in the /—line only accumulates at larger radii where 
the density is lower. The resultant f /i ratio is then obtained from the ratio of the summed 
i— line and /—line fluxes. The key parameter in this distributed f /i emission case is the 
radial location of the densest X-ray source that is capable of contributing to the observed 
i— line and /—line fluxes, i.e., the densest X-ray source that is not heavily attenuated by the 
overlying stellar wind. 

The main difference between these approaches is that the distributed f /i will increase 
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faster with radius than the localized f/i case as shown by Leutenegger et ah (2006). The 
density-squared dependence of the hne emissivity is the primary reason for this shift, al- 
though there is also some change owing to the term in the volume integral. An even more 
complex f/i scenario can be envisioned with the inclusion of a radial temperature distribu- 
tion and stellar wind absorption. This approach is beyond the scope of our current paper, 
and we plan to address this issue in a subsequent paper. 

In this paper (as assumed in all our previous papers on this subject) all quoted values of 
the f/i ratio are based on the localized approach, assuming that the He-like ion Tx is always 
equal to its respective Tl (i.e., the temperature where the given ion's X-ray line emissivity is 
at a maximum which is known for every ion) . We know that a given set of He-like fir lines 
must be forming somewhere in the wind where Tx = Tl and what the localized approach 
provides is the most likely location where this occurs. Furthermore, detailed shock modeling 
such as that of Feldmeier (1995), indicate that the localized approximation for treating f/i 
may indeed be appropriate since the dominant X-ray regions are spatially well separated, 
and the largest contribution to the two lines will arise from the highest density shocked 
plasmas that are capable of producing X-ray emission. Regardless, the main advantage of 
this approach is that it provides results that are only dependent of the adopted atomic 
physics and the UV/EUV model atmosphere. This allows us to study our collection of 
stars in a nearly identical way so as to reveal any significant differences associated with 
basic stellar parameters. However, if the He-like emission is indeed distributed continuously 
throughout the wind as discussed by Leutenegger et al. (2006), the resultant Rfir for a given 
set of He-like fir lines represents the lower radial boundary of the distributed X-ray emission 
integration. Hence, for a wind distribution of X-ray sources there can be no emission below 
this lower radial boundary. In general, these lower radial boundaries are expected to be 
lower than our tabulated localized Rfir values. 

The observed f/i ratios and their associated derived Rfir are given in Tables [3] (MEG 
results) and|l](HEG results). The results for O vil are not given in the HEG table since this 
ion is outside the energy band of the HEG. In addition, although the He- like Ar xvil fir lines 
were detected in Cyg 0B2 No. 8a (Waldron et al. 2004), these results are not presented here 
since none of the other stars have sufficient Ar xvil S/N to carry out a meaningful analysis. 
The most notable observation is that for the supergiants, where there is a clear progression 
in the Rfir being large for low energy ions and small for high energy ions. This behavior 
is what we have referred to as the NSHIP. There is also another interesting observation 
for the supergiants from examination of the Mg xi and Si xill Rfir- We see that essentially 
all of the Rfir for Mg xi are located between 3 and 6 i?*, whereas, the Rfir for Si xill are 
located in a narrower range of 1.8 to 2.3 i?*. For both the giants and MS stars, there is a 
similar dependence, but not as noticeable as for the supergiants. 
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An important observational feature that has emerged from analyses of He-like fir lines 
is the correlation between Rfir and their associated X-ray continuum optical depth unity 
radii {Rr=i) as first noticed by Waldron & Cassinelli (2001). A discussion on the calculation 
of Rr=i is given by Waldron et al. (2004). Figure [H] shows a scatter plot of the dependence of 
Rfir on Rr=i for our luminosity groups (we refer to these plot types as scatter plots because 
the data are shown with error bars in both the ordinate and abscissa directions). The 
correspondence between the two radii is evident in the supergiants which was first shown by 
Waldron & Cassinelli (2002) using a small collection of OB stars. There is a key distinction 
present in all luminosity groups in that there are very few data points that lie fully below the 
dashed line which represents the exact one-to-one correspondence, i.e., Rfir = Rr=i- This is 
consistent with the idea that we do not see any observed line emission arising from below the 
optical depth unity radii which is exactly what one would expect from basic radiation transfer 
arguments, i.e., radiation can only escape from those radial locations where the associated 
optical depths are < 1. For the supergiants, the most likely explanation of this correlation 
is that a very large number of X-ray sources are distributed throughout these winds at 
essentially all radii, and the observed emission line characteristics (line strength and wind 
locations, Rfir) are primarily determined by the dominant X-ray sources. These dominant 
sources are those with the largest emission measures (the density-squared dependence of the 
emissivity) which are no longer hindered by significant wind absorption effects, and their 
maxima emissions arise from their associated X-ray continuum optical depth unity radii 
{Rt=i)- For the giants and MS stars, the X-ray emission is seen to be occurring at wind 
radii that are larger than their associated Rr=i- Clearly, radiation can always escape from 
any region that has a small optical depth, but the fact that there is evidence for radiation 
emerging from radii > Rt=i may be an indication that the number of wind distributed X-ray 
sources in these stars are greatly reduced as compared to the supergiant winds. For example, 
consider a very simple case where there is "one" spherical expanding shock wave, then the 
location of the X-ray emission will always be associated with the radial location of the shock 
wave, independent of the location of Rr=i. By considering say, 1 to 5 expanding shock waves 
at different X-ray temperatures, this may be able to explain the observed scatter shown in 
Figure [6] for the giants and MS stars. Although this is a highly unlikely scenario primarily 
because such a structure would predict strong X-ray variability and we are not aware that 
this is the case. A more likely explanation is that the M of the giants and MS stars may 
actually be larger than is traditionally assumed, so that the Rt=i values may be larger. 

As evident in Figure [6] there are four supergiant data points that do indicate emission 
occurring below Rr=i- One of these is the 5 Ori A Ne ix emission, and the other three 
are associated with the He-like ions of Cyg 0B2 No. 8a. A possible explanation of this 
discrepancy for Cyg 0B2 No. 8a may be related to the uncertainty in the mass loss rate as 
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discussed by Waldron et al. (2004). In the original analysis of 6 Ori A by Miller et al. (2002), 
the derived Ne ix f /i indicated a radial upper limit of ~ 4 i?^,, significantly larger than our 
upper limit shown in Table [31 This discrepancy illustrates the importance of including line- 
blending effects in the extraction of individual line fluxes. The Ne ix wavelength region is 
highly contaminated by many lines. Although we believe that our quoted Ne ix results are 
correct, we still do not have an answer as to why the Ne ix Rfir is ^ 1 R^. A possible answer 
may be related to wind-wind interactions since S Ori A is a well known binary system. It 
may be that the majority of the observed Ne ix emission is actually occurring very close to 
the stellar surface of 6 Ori A binary companion, a BO. 5 III star (Miller et al. 2002). 

A controversy has recently arisen with regards to the Si xill f/i ratio for the late O 
supergiant, ( Ori, as discussed by Leutenegger et al. (2006) and Cohen et al. (2006). In the 
original analysis of ( Ori, Waldron & Cassinelli (2001) used the HEG f/i rather than the 
MEG f/i in their analysis primarily due to possible problems in the interpretation of MEG 
f/i ratio. It is well known that the MEG ancillary response file (ARE) has a significant Si K- 
shell edge (produced by the instrument silicon chip) within the MEG resolution limits of the 
Si XIII /-line which is less prominent in the HEG ARE. Our tabulated ( Ori MEG Si xill 
f/i (see Tables 13] and H]) is consistent with the one derived by Leutenegger et al., but both 
of these are larger than the value determined by Oskinova et al. (2006). However, our HEG 
Si XIII //^ is consistent with the MEG Si xill f/i given by Oskinova et al. (2006). Since most 
of the MEG and HEG determined Si xill f/i ratios are consistent in all the other OB stars, 
it is unclear as to whether this ( Ori discrepancy in the Si xill f/i is related to a specific 
problem with the ( Ori MEG ARE (different software versions), a problem associated with 
the extracted ( Ori count spectrum, or maybe some other unknown problem. For example, 
there are noticeable differences in the Si xill fir lines when comparing the MEG+1 and 
MEG-1 dispersed spectra which could be evidence of contamination by an unknown X-ray 
source in either the MEG+1 or MEG-1 dispersed spectrum in the energy vicinity of the 
Si XIII fir lines. Until this disagreement between the MEG and HEG Si xill f/i for ( Ori 
is resolved, we choose to adopt the HEG f/i for ( Ori as being the most appropriate value 
(as originally proposed by Waldron & Cassinelli 2001). 



5.2. X-ray Temperatures derived from the H/He Line Ratio 

Deriving the temperature distribution versus radius of the X-rays sources in OB stars is 
a major goal of this study. Establishing the X-ray temperatures {Tx) in OB stellar winds is a 
crucial aspect in understanding the mechanisms responsible for the observed X-ray emission. 
Prior to the availability of HETGS data, our knowledge of Tx was limited to one or two 
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temperatures as determined from fitting broad band X-ray spectral data. From HETGS 
data we now know that the OB stellar X-ray emission is produced by a large range of Tx 
(~ 2 to 25 MK), and these temperatures can be used to probe the properties of the X-ray 
sources. If we also have detailed spatial information, we can then obtain a radial distribution 
of the stellar wind temperature structure. In general, since each wind shock may also have a 
temperature stratification which is determined by the extent of the post-shock cooling zone 
(e.g., Feldmeier et al. 1997b), our derived radial distribution of Tx should be construed 
as representing the dominant Tx at each associated radius. By knowing Tx we can also 
determine the pre-shock velocity relative to the shock front (Urei), defined as Urei = V — Vs 
where V is the pre-shock gas velocity, Vs is the velocity of the shock front, and V and Vs are 
measured in the rest frame of the star. From the Rankine-Hugoniot relation, the dependence 
between Tx (or post-shock temperature) and Urei is given by 



The quantity Urei is the fundamental parameter required to determine the physical 
characteristics of a shock. In general, Urei can be used to determine the shock strength, 
post-shock temperature, and relative post-shock speed ( = 1/4 Urei)- For a basic shock 
model description, since the pre-shock gas velocity in the rest frame {V) is defined as equal 
to the ambient wind velocity (hereafter denoted as Vq), Urei can also be used to determine 
the speed of the shock front and a constraint on the shock front radial location provided we 
know the radial dependence of Tx- It is implied that all of these quantities ( Tx, Urei, V, 
Vs, and Vq) are dependent on radius. 

The relevance of the H/He emission line ratio as a diagnostic of the X-ray source 
temperatures in OB stellar winds has been explored by Schulz et al. (2000), Miller et al. 
(2002), and Waldron et al. (2004). The advantages of using H/He ratios are: 1) the H- 
like and He-like lines are very strong in HETGS data; 2) the line ratios display a strong 
dependence on the X-ray temperature as shown in Figure [5l and; 3) line-blending effects on 
the H/He ratio are minimal, with the exception being the Neon H/He ratio. We point out 
that our results shown in Figure [5] are slightly different from those shown by Miller et al. 
(2002) in that their He-like line emission used in their H/He ratio included all the fir lines, 
while we choose to only consider the He- like r— line to represent the He-like line emission 
(see discussion by Waldron et al. 2004). 

The one possible disadvantage in the usage of the H/ He ratio with regard to studies of 
OB stars is that the H-like and He-line ions could possibly form in different regions of the 
stellar wind which implies that they may suffer from different amounts of wind absorption 
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and/or resonance line scattering. For example, Porquet et al. (2001) suggest that the 
observed G-ratio may be larger than expected since the He-like r— line can be strongly 
affected by resonance hne scattering escape probability effects in contrast to the i— and 
/—lines. Hence, the presence of resonance line scattering would lead to an underestimate 
of the actual X-ray temperature. In the Appendix we explore the possible effects of X- 
ray continuum absorption and resonance line scattering on the observed H/He ratios. Our 
results suggest that the temperatures derived from H/He ratios are good indicators of the 
true X-ray temperatures. Although the expected temperature dependent behavior of this 
ratio shown in Figure [5] is based on the MEG energy resolution limits to determine the 
range of line-blending, we find that the HEG energy resolution produces only very minor 
differences. 

We extract H/He temperatures, Thhe, for all available OB H-He line pairs using the 
expected temperature dependent H/ He ratios shown in Figure [5] to determine the range in 
temperature associated with the range in the observed ratios. Prior to the extraction of 
Thhb, all lines in these observed ratios are corrected for ISM absorption, but no attempt 
is made to estimate wind absorption and resonance line scattering. The observed H/He 
ratios and their inferred Thhs are listed in Tables O^MEG results) and [H] (HEG results). For 
completeness, the observed MEG and HEG G-ratios are listed in Table [71 

The Thhe data listed in Tables [5] and [6] illustrate several important results: 1) these X- 
ray sources show a large range in X-ray temperature, from ~ 2 to 23 MK; 2) all of the Thhe 
are found to be within the temperature range as specified by their respective H and He T^, 
i.e., the Tl range shown in Tables [5] and [6] where Tl represents the temperature associated 
with the maximum line emission of a particular ion (see Table [2]); 3) the consistency between 
Thhc and for all ions can be interpreted as a verification that the OB observed X-ray 
emission lines arise from a thermal plasma since Ti represents the collisional ionization 
equilibrium temperature associated with a given line's maximum line emission; 4) for a 
given ion and a given luminosity class, the values of the derived Thhe are very similar, and; 
5) comparisons of the mean Thhb for each luminosity group show no significant differences, 
except possibility the supergiant Mg xi mean Thhe which is somewhat larger than the other 
luminosity groups, but this difference is traced to the higher Thhs of the Cyg 0B2 stars. 
Note that the mean Thhe for the MS stars does not include the results for 6^ Ori C, a well 
known peculiar magnetic star. The most obvious difference between the 6^ Ori C Thhc for 
each H/ He line pair is that they are larger than those of the other stars in our sample (see 
Tables & E]) • In addition, the 6^ Ori C temperatures are only slightly larger than those of 
T Sco which has also been confirmed to have magnetic structures (Donati et al. 2006). 

If the H/He line pairs were sensitive to wind attenuation owing to continuum opacity. 
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we would have expected a wider dispersion in Thhc since our program stars contain a wide 
diversity in stellar wind properties. However, there are a few stars in each luminosity group 
where there is some noticeable difference among individual star Thh^. with respect to the 
mean of the group. This may signify certain effects such as, an overall higher tempera- 
ture structure, wind absorption, and resonance line scattering. In general, since the H/He 
temperature diagnostic appears to predict temperatures that appear to be independent of 
wind structure (i.e., individual stars), we suggest that this ratio should be considered to 
be the best available diagnostic for establishing the X-ray source region temperatures. The 
overall consistency in these temperatures is a very interesting result, in that, regardless of 
luminosity group, as well as individual stars, the expected temperature for any given H-like 
to He-like line pair is always the same which implies a global commonality among the X-ray 
temperature distributions for all OB stars. 



6. The Stellar Wind Distributions of HWHM, Vp, Tx, and Urei 

In the previous sections, we have determined the line emission parameters {HW H M , Vp, 
and EMx), spatial locations of the He-like ions (Rfir), and X-ray temperatures {THHe)- We 
now combine this information to find the stellar wind spatial distribution of the line emission 
parameters and temperatures. We will not explore the radial dependence of EMx primarily 
because these are the observed emission measures and not the intrinsic emission measures, 
plus there is some uncertainty owing to the makeup of EMx such as shock density and 
temperature, and geometric factors, such as surface area and thickness. In our discussions 
of the radial dependence of HWHM, Vp, and Urei it is beneficial to use "velocity versus 
velocity" type plots where the abscissa for any given radius, R, is defined as Vo(-R)/foo- 
The transformation between R and Vq is obtained through the standard "/3-law" velocity, 
Vo{R)/voo = (1 - Ro/Rf, assuming a /3 = 0.8. 

Figure [7] shows scatter plots of the normalized HWHM and Vp of the He-like ions versus 
the normalized ambient wind velocities [Vo(-R/jr)]; where all quantities are evaluated at their 
Rfir for all OB luminosity classes. From the basic shock model description one would expect 
the HWHM and Vp at any given radius (R) to be less than the Vo{R) associated with the 
line formation region since there is a significant decrease in the flow velocity across a shock. 
For lines forming at large radii, the HWHM are < VoiRfir) as shown in the upper panel of 
Figure [3 However, for lines near the star, although the HWHM are < f oo these HWHM 
are > Vq at these low radial positions, which presents a problem in that the lines formed near 
the star are too broad. This is also a problem for the de-shadowing shock model (Feldmeier 
et al. 1997a) which predicts that the post-shock velocities are comparable to the ambient 
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wind velocities. Our results indicate that the spread in HWHM ranges from to 0.6 foo at 
all radii. In the lower panel of Figure [71 we see that the Vp values are barely shifted and the 
spread is only ±0.2 Voo- Therefore, we find that observed distributions of both the HWHM 
and Vp are nearly independent of the wind location of the X-ray emitting plasma. This is 
an odd result that has only become apparent from our multi-star study, i.e., there was no 
indication of such behavior emerging from analyses of individual stars. 

Figure [H] shows a scatter plot of the dependence of Thhe on Rjir for the three luminosity 
groups. Our results confirm the initial results reported by Waldron (2005), the Thhe of the 
supergiants show a well defined radial distribution that decreases outward from the stellar 
surface as evident in both the MEG and HEG data. The supergiant data show a strong 
correlation in which the highest temperatures (20 MK) only occur very close to the star, and 
the lowest temperatures (2 MK) occur only in the outer wind regions. Although this behavior 
is not as obvious in the giants and MS stars, it is clear that the highest temperatures are only 
located near the surface, but low temperature can occur anywhere within the wind. It is 
also interesting that the observed radial distribution of Thhc indicates that for all luminosity 
classes there is a radial dependent maximum X-ray temperature which decreases with radius. 
To our knowledge, this behavior has not been predicted by any shock model. The most 
likely reason as to why the supergiants show such a tight correlation between temperature 
and radius is that the supergiants have larger wind densities and column densities, hence, 
we cannot see the low temperature regions that are present at small wind radii. This follows 
from the optical depth unity argument. For the giants and MS stars shown in Figure [HI 
one sees a large range in temperature at small and intermediate radii with no cutoff as seen 
in the supergiant case. This is probably also the case for the supergiants as well, but we 
cannot see the lower temperatures at small and intermediate radii due to larger stellar wind 
absorption. Furthermore, since all three luminosity groups show high energy ions existing 
very close to the stellar surfaces, the NSHIP is a common feature identified with all OB 
stars. 

We now choose to examine the radial dependent distribution of the He-like ion values of 
Urei (the relative pre-shock velocity discussed in Sec. 5.2) with regards to the radial depen- 
dent ambient wind velocity, Vq, where Urei and Vq are evaluated at their Rfir- The values 
of Urei are obtained directly from eq. \5\ using the radial dependent X-ray temperature struc- 
tures [THHe{Rfir)] showu in Figure [HI Our discussion will incorporate a useful parameter, r], 
which in general should also be dependent on radius, and is defined as 
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This ratio measures the ratio of the pre-shock velocity relative to the shock front to the 
fixed frame pre-shock velocity. For outflows {V > and Vs > 0) rj has a maximum value of 
1 when Vs = which implies that the shock front is stationary in the rest frame of the star, 
and Urei = V, the maximum value of Urei- For the case of in-falling gas clumps, i.e., Vs < 0, 
1] can be > 1, and the actual value of rj is determined by the magnitude of the clump in-fall 
velocity. 

In the following discussion we will address the results as relevant to the basic shock model 
description for an outflowing gas, then, by definition, V = Vq (the ambient wind velocity). 
The scatter plot of Urei versus Vq evaluated at the associated MEG determined Rfir is shown 
in Figure [9] for all the OB luminosity classes. The data display several interesting features: 1) 
there does not appear to be any obvious correlation between Urei and Vq; 2) the majority of 
the data indicate a range in Urei/voo from ^ 0.2 to 0.5; 3) the majority of the lines forming at 
large radii (i.e., with Vo{Rfir)/voo > 0.6) occur where r] < 0.5, indicating weak relative pre- 
shock velocities which implies that the shock front velocity (Vs) is becoming comparable to 
the pre-shock velocity (V^), since = V{1 — rj), as r] 0, then, Vs ^ V; 4) at intermediate 
ambient wind speeds (~ 0.5 Voo, a radius of ~ 1.7 /?=„), < 1 which means that Vs ^ 0, 
i.e., the shock front is stationary in the star's rest frame, and; 5) the most interesting result 
is associated with those lines that form very near to the stellar surface where Vq < 0.2 v^o, 
or a radius < 1.15 i?*, with a large range in Urei from 0.2 to 0.8 Voo, and these are the lines 
discussed earlier regarding the near-star-high- ion problem (NSHIP). 

The data appear to support the idea that the majority of the He-like ion lines are 
consistent with a basic shock model interpretation provided that these winds have shock front 
velocities that are initially stationary at various points in the wind where t] = 1 and then they 
accelerate outwards where eventually the shock front velocities become comparable to the 
pre-shock velocities and these weak shocks are no longer capable of producing detectable X- 
ray emission. This interpretation is consistent with our derived X-ray temperature structures 
shown in Figure [H It is also likely that similar conclusions may be attainable from other 
shock models with appropriate parameter adjustments. 

However, the He-like ion lines at low radial locations {Vo{Rfir) < 0.2 Voo) cannot be 
explained by this basic shock model description. We explore two possible alternatives, and 
use the observed data point with Urei = 0.5 and Vq = 0.1 as our sample for comparisons. For 
an in-falling clump, we find an rj = 5 and an in-fall velocity of Vs = —0.4 Voo- It is unclear 
as to whether such large in-fall velocities can occur at small radii. Now we examine whether 
the rapid acceleration in the de-shadowing instability model is applicable. As discussed 
earlier, for outflowing gas, the maximum of 77 is 1 which means that V = Urei = 0.5foo, 
and this acceleration must occur at a radius < 1.1 i?*. From inspection of several numerical 
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simulations (Owocki et al. 1988, MacFarlane & Cassinelli, 1989, Feldmeier 1995, Feldmeier 
et al. 1997a, Runacres & Owocki 2002) we do not see any evidence for large pre-shock 
velocities of O.Svoo at such low stellar radii. 



7. Summary 

The main goal of this paper has been to present a large collection of observational 
results and examine these in a uniform way so as to draw out general facts regarding OB 

stellar X-ray emission characteristics. We have analyzed the Chandra HETGS MEG and 
HEG data from 17 OB stars. Although the HEG spectra are typically much weaker than 
their associated MEG spectra, we have chosen to use both data sets to illustrate the rather 
good agreement between their observed and derived attributes. As a slight departure from 
our original goal of analyzing only "normal OB stars" we have also included Ori C in our 
sample, primarily to allow one to compare the X-ray properties of a know peculiar 0-star 
with those of other OB stars. 

We use a well recognized Gaussian line-fitting method to extract the pertinent X-ray 
emission line parameters which ensures that our results are easily reproducible. This model- 
independent line fitting procedure has allowed us to determine the observed line emission 
parameters [HWHM, Vp, EMx, and line fiux) and fine emission ratios {f/i and H/He) in 
order to obtain an easily verifiable description of these X-ray emission properties. We have 
primarily focused on searching for luminosity class regularities in the data, and the radial 
distributions of several X-ray derived parameters. 

HWHM/voo Histograms: The HWHM /v^ plots show peaks that are well below the 
wind terminal velocity for all luminosity classes. We have argued that the differences from one 
luminosity class to another in regards to the asymmetries in these distributions are due to the 
differences in wind column densities. The fact that a large percentage of HWHM /voo < 0.5 
indicates that the majority of the observed line emissions occur relatively deep in the wind 
where the wind flow is accelerating. 

Vp/voo Histograms: The Vp/voo plots show that the majority of the lines are symmetric, 
with very little line-shifts. There is some tendency for stars of all luminosity classes to show 
a small but finite blue-ward asymmetry in their line shift distributions. Such an asymmetry 
would occur if the X-rays are primarily arising from the near side of the star. Several ideas 
have been discussed, ranging from in-falling clumps, the orientation of elongated clumps 
with respect to observer's LOS (line-of-sight), a two-component wind geometry consisting 
of a polar wind and an equatorial wind which have different wind densities and velocity 
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structures, to a major decrease in mass loss rates as compared to traditional values. The 
actual explanation is not yet clear. 

Emission Measure Histograms: The supergiant EMx values show a well defined peak at 
log EMx = 54.75 which is due to the fact that these stars all have similar M, Voo, and R^^. 
However, this is not the case for the giants and MS stars, and the resultant distributions are 
more spread out. In addition, the MS stars show a double peaked distribution in EMx which 
is due to the fact that there is a broader range in wind properties along the MS spectral 
class. 

X-ray Source Locations: Using the He-like f /i line ratios we have derived the stellar wind 
location, Rfir, of the O, Ne, Mg, Si, and S He-like ion emission. The results clearly emphasize 
the NSHIP (near-star-high-ion problem) as evident from the S XV derived Rfir- We stress 
that these Rfir are based on the localized interpretation and a distributed interpretation 
would result in a distribution that starts at an even smaller Rfir- Our result confirms what 
has been shown for individual stars, all OB stars display the same basic distributions with 
the He-like O and Ne ions located in the outer wind regions, the He-like Mg and Si ions 
at intermediate locations, and the He-like S ions located near the stellar surface. For all 
supergiants, the MEG He-like Mg and Si f/i ratios predict a narrow radial region of 3.7 to 
5.6 i?*, and 1.5 to 3.0 R* respectively, whereas the He-like Ne and O f /i ratios predict radial 
locations from 6.2 to 13.5 -R*. The MEG He-like S f/i predict locations that are essentially 
on the surface. In general, there is good agreement between the MEG and HEG derived 
Rfir, with one exception, the Si xill results for ( Ori (see discussion in Sec. 5.1). 

X-ray Temperatures: We think our most interesting results are in regards to the Tx 
radial distributions. We have derived X-ray temperatures (THHe) using the H/He line ratio 
temperature diagnostic. First, we emphasize the Thhc should be considered as an average 
temperature of the H-like and He-like ions (i.e., basically an average of the associated ion 
peak Tl given in Table [2]). Overall our results are consistent with this interpretation, e.g., 
by comparing the Tl range with the mean values listed in Tables [5] and [61 In general, for 
a given H/He line pair, the resultant Thhe are found to be essential the same regardless 
of the luminosity class. These results can provide valuable information regarding the shock 
formation processes that are actually operating in OB stellar winds. However, there are 
some notable differences. The silicon mean Thhe for all luminosity classes are at the lower 
end of the Tl range with a larger spread in Thhc as compared to the other H-He line pairs. 
The spread in Thhc suggests that the mechanism producing these higher temperatures may 
be dependent on the individual stellar characteristics and probably related to resonance line 
scattering as discussed in the Appendix. The three stars that have their silicon Thhe closer 
to the upper Tl limit are the two Cygnus 0B2 stars and the known peculiar star, 9^ Ori C. 
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As to whether these Cygnus 0B2 stars are similar to 9^ Ori C remains to be determined. 
Ideally, one would like to have a more complete collection of sulfur T^He to explore the higher 
temperature behavior, but the data do not allow reasonable signal-to-noise extractions. 

The Correlation Between Rfir and Rt=\: Since the initial detailed studies of O super- 
giants (e.g., Waldron & CassineUi 2001; CassineUi et al. 2001), the basic idea that the Rfir 
arc correlated with Rr=i has continually surfaced as a reasonable observational result. To 
test this correlation we have examined the luminosity class dependence. Our results show 
that this correlation is clearly evident in the supergiants. But, the data from the giants and 
MS stars indicate that the X-ray emission arises primarily from radial locations above their 
Rr=i surfaces. We suggest two possible explanations, either the giants and MS stars have 
larger mass loss rates, or their number of distributed wind shocks are significantly reduced as 
compared to those in supergiant winds (see discussion in Sec. 5.1). However, what is clearly 
emerging from our results is that for all luminosity classes, there are very few observed X-ray 
sources where Rfiris less than Rt=i. This supports our basic interpretation of this correla- 
tion. One can only see those X-rays that are capable of escaping the wind, and since the line 
emission is proportional to rig^, one is likely seeing the peak of the contribution function as 
deep in the wind as possible. If the mass loss rates of these stars are indeed much smaller 
than previously thought, one would have expected to see an uncorrelated scatter plot, i.e., 
we should also be seeing a scattering of X-ray emission from regions where Rfir < Rt=i. We 
do not see any plausible explanation as to why such a correlation should hold in a highly 
porous or clumped wind since the X-ray continuum absorption is determined by the "cool" 
stellar wind opacity (as determined by bound-free transitions) which is linearly dependent 
on the density in the X-ray energy range. Hence, a reduced wind density or clumped wind 
structure would predict Rr=i values much lower than those adopted in our study. In addi- 
tion, this correlation implies that we can now predict a particular He-like f/i ratio by simply 
calculating the expected Rt=i- We therefore choose to present this as a challenge to the idea 
that the M values are so much lower, and perhaps also a problem for a clumped wind. 

Radial Distribution of HWHM and Vp.' Since we know the radial locations of the 
He-like fir emitting ions we can explore the radial dependence of the He-like X-ray line 
properties. Almost all line HWHM are found to be less than the radial ambient wind 
velocity as expected from all shock model descriptions, but there is no apparent correlation 
between HWHM and the ambient wind velocity. However, several deeply embedded sources 
have HWHM that are significantly greater than the radial ambient wind velocity which 
implies a new interesting problem, the lines from deep in the flow are too broad! There is 
no obvious radial dependence of Vp where the majority of observed Vp/voo remain between 
±0.2 regardless of wind location, a totally unexpected result. This radial dependence of Vp 
is an important problem that must be addressed before we can fully understand the source 
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of the X-ray emission from OB stars. 

Radial Distribution of X-ray Temperatures: Although in principle we can only address 
the radial locations of temperatures associated with the He-like fir emission lines, we have 
argued that the H/ He line ratios derived Thhs provide reasonable estimates of the average 
temperature associated with the H-like and He-hke ions. Waldron (2005) was the first to 
show a correlation between the Thhc and Rfij. for OB supergiants. We have now shown 
that this correlation is consistent with both the MEG and HEG results. Our analysis has 
verified that the temperature steadily decreases outward through the wind with the highest 
temperatures occurring near the stellar surface. For the other luminosity classes, although 
the tight correlation breaks down, we still find that the highest temperatures only occur 
near the stellar surface. Hence, there is no evidence of any high temperature in the outer 
wind regions which is perhaps surprising since shock velocity jumps could in principle be 
at their largest at these large radial locations. Our results support two interesting results. 

1) There appears to be a well defined radial dependent maximum X-ray temperature. This 
radial dependent maximum temperature can be extremely useful in determining basic shock 
characteristics such as, the efficiency of the conversion of shock energy into X-ray emission. 

2) The magnitude of the derived X-ray temperatures and associated radial distributions are 
dependent on wind density. In the dense winds of the supergiants we see little evidence of 
low Tx at low radial locations, whereas, in the lower density winds of the giants and MS 
stars, a large range in Tx exists at lower radii. The implications are clear. X-ray emission 
is probably existent throughout these winds at all temperatures provided they are below 
the radial dependent maximum temperature. However, as the wind density gets larger, one 
cannot detect the low Tx at low radii due to wind absorption effects. This is exactly what 
one would expect based on our optical depth unity radii arguments. Furthermore, for all OB 
luminosity classes we find that the radial dependence of the pre-shock velocities relative to 
the shock front {U,.ei) evaluated at T^He, indicate that the majority of all lines at large radial 
distances are consistent with the basic shock model interpretation. These results also support 
outward accelerating shock front velocities which are consistent with the decreasing X-ray 
temperature structure. However, a fundamental problem exists at low radial locations (< 1.2 
i?*) where there are several lines, characterized by the high ion stages, with temperatures 
that appear to be too hot, and line profiles that have HWHM that are too broad. As to 
whether any wind shock mechanism can explain these too hot, and too broad X-ray emission 
lines, needs to be investigated. 

Some of the problems discussed here have recently been addressed in repeated analyses of 
the X-ray emission from a few select 0-stars such as C, Pup and C, Ori using a semi-empirical 
shock model description (e.g., Oskinova et al. 2006; Leutenegger et al. 2006; Cohen et 
al. 2006). Cohen et al. (2006) argues that all emission lines from C, Ori can be explained 
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by a shock model description if the X-rays are distributed above 1.5 i?*, provided that the 
wind density is greatly reduced or highly clumped. Cohen et al. further argue that there 
is no problem with understanding the X-ray emission from OB stars. However, we see two 
potential problems with these conclusions: 1) their modeling efforts assume that everywhere 
Tx = Tl regardless of the ion being considered, and as we have shown in Section 5.2, this 
is not a good assumption, and more importantly; 2) the X-ray emission from below 2 i?* 
is a notoriously difficult problem due to the line drag effect (Lucy 1984) or strong source 
function gradients (Owocki & Puis 1999) that affect the hne-driven instability. Hence, we 
argue that until we have a better understanding of the shock properties below 2 i?*, such 
conclusions are premature. 

We have provided a summary of our results and emphasized the critical problems asso- 
ciated with current interpretations. Our results have verified that essentially all OB X-ray 
emission lines are un-shifted. We have introduce a new problem associated with the X-ray 

temperatures and spatial locations which we have labeled as the "near-star high-ion problem " 
{NSHIP) and we believe this problem is far more critical than the un-shifted line prob- 
lem. Resolution of this problem will lead to a better understanding of the X-ray production 
mechanisms at work in OB stars. It seems that the wind shock model may be responsible 
for the observed X-ray emission in the outer wind regions (although we still must figure out 
why these lines are un-shifted), but we must consider alternatives for the remaining X-ray 
emission. The un-shifted line problem has been addressed by either lowering the mass loss 
rates or proposing a clumpy wind. Either approach produces the same desired effect, an 
overall reduction in the stellar wind X-ray absorption. This reduction in X-ray absorption 
allows for more X-rays to emerge from the far side of the star (the red-ward emission) which 
makes it possible to explain symmetrical line profiles. Although this provides an explanation 
for the observed line symmetry, we arc then left with a dilemma as to why such a strong 
correlation exist between Rfi,. and Rr=i (as determined from the traditional mass loss rates) 
where the premise of this correlation is based on fundamental radiation transfer arguments. 

With regards to the observed X-ray emission not explainable by standard wind shocks, 
a continuing possibility is that these deeply embedded X-ray sources are associated with 
magnetic fields on or near the stellar surface. It is becoming increasing clear that these stars 
are theoretically allowed to have surface magnetic structures (e.g., MacGregor & CassineUi 
2003; MuUan & MacDonald 2005). In any case, if there is an alternative source of X- 
rays operating in these stars, it will be an exciting new feature associated with OB stellar 
astronomy. 

We would like to thank the anonymous referee for his detailed critique of our manuscript 
and several informative suggestions. WLW acknowledges support by award GO2-3027A 
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TM3-4001 issued by the Chandra X-ray Observatory Center. Chandra is operated by the 
Smithsonian Astrophysical Observatory under NASA contract NAS8-03060. 

A. Wind Absorption and Resonance Line Scattering Effects 

To examine the possible effects of continuum absorption and escape probability effects, 
consider the following simple illustration. If H represents the observed H/ He ratio and Hq 
represents the intrinsic temperature dependent H/He ratio, then, 

H = Ho \ "^^'^ exp ( THe - th ) (Al) 
1 + apH 

where thc and th are respectively the X-ray continuum wind absorption optical depths for 
the He-like and H-like lines, and pne and pn represent the associated resonance line scattering 
T using the approximate escape probability formalism developed by Osterbrock (1974) (a is 
a constant equal to 0.58). First we examine the effects of continuum absorption only (i.e., 
PHe = Ph = 0). Since the OB stellar wind opacities scale roughly as A'^, all H-He line pairs 
observable in HETGS data are expected to have thc > th if both lines are formed at the 
same radial position. Hence, for the case of negligible resonance line scattering, H would 
be expected to be larger than Hq which means that a temperature derived from H will 
be greater than the actual value. There is one exception to this scaling. For the oxygen 
H-He lines, the reverse is expected, i.e., th > T^e which would imply that the H derived 
temperature is actually lower than the real value (from Table [2] notice that the aw of O vill 
is significantly larger than aw of O vil). Now we consider the presence of finite resonance 
line scattering. If we again assume that both H-like and He-like lines are formed at the same 
density and temperature, then ph^ = 2pH based on atomic physics considerations. Hence, 
an additional increase in H would also occur for the case of finite resonance line scattering. 
For example, consider an extreme case where ths - th = ^ and pn = 1 {PHe = 2), then H 
= 3.7 Hq- Now suppose we know that are observing a Si plasma at a temperature of 10 
MK [Hq = 0.4), then under these conditions the observed if = 1.5 (FiglS]). The resultant 
observed H-He ratio would predict a temperature ~ 15 MK respectively. However, a more 
likely scenario would be thc - th = since it seems that all observed line emissions seem to 
always arise from near their respective X-ray continuum optical depth unity radii as shown in 
Figure O For this case, the H would be equal to 1.4 Hq and the predicted H-He temperature 
would be ~ 11 MK, much closer to the actual temperature. Hence, we argue that these H/ He 
derived temperatures provide very good estimates of the actual X-ray temperatures. 



-38- 



REFERENCES 

Berghdfer, T. W., Schmitt, J. H. M. M., & Cassinelli, J. P. 1996, A&ASS, 118, 481 
Berghofer, T. W., Schmitt, J. H. M. M., Banner, R., & Cassinelli, J. P. 1997, A&A, 322, 167 
Blumenthal, G. R., Drake, G. W. F., & Tucker, W. H. 1972, ApJ, 172, 205 
Cargill, P. J., & Pneuman, G. W. 1984, ApJ, 276, 369 
Cash, W. 1979, ApJ, 228, 939 

Cassinelh, J. P., Cohen, D. H., MacFarlane, J. J., Sanders, W. T., & Welsh, B. Y. 1994, 
ApJ, 421, 705 

Cassinelh, J. P., & Olson, G. L. 1979, ApJ, 229, 403 

Cassinelli, J. P., Waldron, W. L., Sanders, W., Harnden, F., Rosner, R., & Vaiana, G. 1981, 
ApJ 250,677 

CassineUi, J. P., Miller, N. A., Waldron, W. L., MacFarlane, J. J., & Cohen, D. H. 2001, 
ApJ 554, L55 

CassineUi, J. P., & Swank, J. H. 1983, ApJ, 271, 681 

Cohen, D. H., Cooper, R. G., MacFarlane, J. J., Owocki, S. P., Cassinelli, J. P., & Wang, P. 
1996, ApJ, 460, 506 

Cohen, D. H., de Messiees, G. E., MacFarlane, J. J., Miller, N. A., Cassinelli, J. P., Owocki, 
S. P., & Liedahl, D. A. 2003, ApJ, 586, 495 

Cohen, D. H., Leutenegger, M. A., Grizzard, K. T., Reed, C. L., Kramer, R. H., & Owocki, 
S. P. 2006, MNRAS, 368, 1905 

Donati, J. -F., Howarth, I. D., Jardine, M. M., Petit, P., Catala, C, Landstreet, J. D., 
Bouret, J. -C, Alccian, E., Barnes, J. R., Forveille, T., Paletou, F., & Manset, N. 
2006, MNRAS, 370, 629 

Feldmeier, A. 1995, A&A, 299, 523 

Feldmeier, A., Kudritzki, R. -P., Palsa, R., Pauldrach, A. W. A., & Puis, J. 1997b, A&A, 
320, 899 

Feldmeier, A., Puis, J., & Pauldrach, A. W. A. 1997a, A&A, 322, 878 
Feldmeier, A., Shlosman, I., & Hamann, W. -R. 2002, ApJ, 566, 392 
Feldmeier, A., Oskinova, L., & Hamann, W. -R. 2003, A&A, 403, 217 
Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1039 
Gabriel, A. H., & Jordan, C. 1969, MNRAS, 145, 241 



-39- 



Gagne, M., Oksala, M. E., Cohen, D. H., Tonnesen, S. K., ud-Doula, A., Owocki, S. P., 
Townsend, R. H. D., & MacFarlane, J. J. 2005, ApJ, 628, 986 

Groenewegen, M. A. T., Lamers, H. J. G. L. M., & Pauldrach, A. W. A. 1989, A&A, 221, 
78 

Harnden, F. R., Jr., Branduardi, G., Gorenstein, R, Grindlay, J., Rosner, R., Topka, K., 
Elvis, M., Pye, J. P., & Vaiana, G. S. 1979, ApJ, 234, L51 

Howarth, I. D., Bolton, C. T., Crowe, R. A., Ebbets, D. C, Fieldus, M. S., FuUerton, A. 
W., Gies, D. R., McDavid, D., Prinja, R. K., Reid, A. H. N., Shore, S. N., & Smith, 
K. C. 1993, ApJ, 417, 338 

Howarth, I. D., & Prinja, R. K. 1989, ApJS, 69, 527 

Howk, J. C, CassineUi, J. P., Bjorkman, J. E., & Lamers, H. J. G. L. M. 2000, ApJ, 534, 
348 

Hubeny, L, & Lanz, T. 1995, ApJ, 439, 875 
Ignace, R., & Gayley, K. 2002, ApJ, 568, 594 

Kahn, S. M., Leutenegger, M. A., Cottam, J., Rauw, G., Vreux, J.-M., den Boggende, A. J. 
F., Mewe, R., & Giidel, M. 2001, A&A, 365, L312 

Koch, R. H., & Hrivnak, B. J. 1981, ApJ, 248, 249 

Kramer, R. H., Cohen, D. H., & Owocki, S. P. 2003, ApJ, 592, 532. 

Kurucz, R. L. 1993, in ASP Conf. Ser. 44, Peculiar versus Normal Phenomena in A- 
Type and Related Stars, eds., M. M. Dworetsky, F. CasteUi, & R. Faragglana (San 
Francisco: ASP), p. 87 

Lamers, H. J. G. L. M., & Leitherer, C. 1993, ApJ, 412, 771 

Lamers, H. J. G. L. M. & Morton, D. C. 1976, ApJS, 32, 715 

Leutenegger, M. A., Paerels, F. B. S., Kahn, S. M., & Cohen, D. H. 2006, ApJ, 650, 1096 

Leutenegger, M. A., Owocki, S. P., Kahn, S. M., & Paerels, F. B. S. 2007, ApJ, in press 

Leitherer C. 1988, ApJ, 326, 356 

Lucy, L. B. 1984, ApJ, 284, 351 

Lucy, L. B. 1982, ApJ, 255, 286 

Lucy, L. B., & Solomon, P. M. 1970, ApJ, 159, 879 

Lucy, L. B., & White, R. L. 1980, ApJ, 241, 300 

MacFarlane, J. J., & CassineUi, J. P.1989, ApJ, 347, 1090 



-40- 



MacFarlane, J. J., Cassinelli, J. P., Welsh, B. Y., Vedder, P. W., Vallerga, J. V., & Waldron, 
W. L. 1991, ApJ, 380, 564 

MacGregor, K. B., & Cassinelli, J. P. 2003, ApJ, 586, 480 

Mewe, R., Raassen, A. J. J., Cassinelli, J. P., van der Hucht, K. A., Miller, N. A., & Gudel, 
M. 2003, AdSpR, 32, 1167 

Miller, N. A., Cassinelli, J. P., Waldron, W. L., MacFarlane, J. J., & Cohen, D. H. 2002, 
ApJ, 577, 951 

MuUan, D. J., & MacDonald, J. 2005, MNRAS, 356, 1139 
MuUan, D. J., & Waldron, W. L. 2006, ApJ, 637, 506 

Ness, J. -U., Brickhouse, N. S., Drake, J. J., & Huenemoerder, D. P. 2003, ApJ, 598, 1277 

Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2006, MNRAS, 372, 313 

Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2004, A&A, 422, 675 

Osterbrock, D. E. 1974, Astrophysics of Gaseous Nebulae, (San Francisco, CA: W. H. Free- 
man & Company), p. 83 

Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 

Owocki, S. P., & Cohen, D. H. 2006 ApJ, 648, 565 

Owocki, S. P., & Cohen, D. H. 2001 ApJ, 559, 1108 

Owocki, S. P., & Puis, J. 1999 ApJ, 510, 355 

Porquet, D., Mewe, R., Dubau, J., Raassen, A. J. J., & Kaastra, J. S. 2001, A&A, 376, 1113 
Runacres, M. C, & Owocki, S. P. 2002, A&A, 381, 1015 
Rybicki, G. B., & Hummer, D. G. 1978, ApJ, 219, 654 

Savage, B. D., Bohlin, R. C, Drake, J. F., & Budich, W. 1977, ApJ, 216, 291 

Schulz, N. S., Canizares, C. R., Huenemoerder, D., & Lee, J. C. 2000, ApJ, 545, L135 

Schulz, N. S., Canizares, C. R., Huenemoerder, D., & Tibbcts, K. 2003, ApJ, 595, 365 

Seward, F. D., Forman, W. R., Giacconi. R., Griffiths, R. E., Harnden, F. R., Jr., Jones, C, 
& Pye, J. P. 1979, ApJ, 234, L55 

ShuU, J. M., & Van Steenberg, M. E. 1985, ApJ, 294, 599 

Smith, R. K., & Brickhouse, N. S. 2000, Rev. Mexicana Astron. Astrofis. Ser. Conf., 9, 134 
Smith, R. K., & Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJ, 556, L91 
ud-Doula, A., & Owocki, S. P. 2002, ApJ, 576, 413 



-41 - 



Waldron, W. L. 1984, ApJ, 282, 256 

Waldron, W. L. 2005, AIP Conf. Ser., "ITAMP Workshop: X-ray Diagnostics of Astrophys- 
ical Plasmas", (held in Cambridge, MA, Nov 2004), 774, p. 353 

Waldron, W. L., & CassineUi, J. P. 2001 ApJ, 548, L45 

Waldron, W. L., & CassineUi, J. P. 2002, ASP Conf. Ser., The High Energy Universe at 
Sharp Focus: Chandra Science, Vol. 262, eds., E. M. Schlegel & S. Dil Vrtilek, (San 
Francisco), p. 69 

Waldron, W. L., CassineUi, J. P., Miller, N. A., MacFarlane, J. J., & Reiter, J. C. 2004, ApJ, 
616, 542 

Waldron, W. L., Corcoran, M. F., Drake, S. A., & Smale, A. P. 1998, ApJS, 118, 217 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-42 - 



Table 1. Adopted Stellar Parameters and HETGS Data Identification Numbers 



Star 


Spectral 


HETGS 


d 




R* 


M 




Nwo 




Type 


Obs ID 


kpc 


K 


Rq 


lO-f^Mgyr-l 




1022cm-2 


Supergiants (I, II) 


















CPup 


04 If 


640 


0.43 


42400 


16.5 


2.40 


2200 


1.98 


Cyg OB2 No. 9 


05 f 


2572 


1.82 


44700 


34.0 


12.70 


2200 


5.10 


Cyg OB2 No. 8a 


05.5 1(f) 


2572 


1.82 


38500 


27.9 


13.50 


2650 


5.48 


<5 Ori A 


09.5 II 


639 


0.50 


32900 


17.0 


1.07 


2300 


0.82 


C Ori A 


09.7 lb 


610, 1524 


0.50 


30900 


31.0 


2.50 


2100 


1.15 


e Ori 


BO la 


3753 


0.46 


28000 


33.7 


4.07 


1500 


2.42 


Giants (IV, III) 


















HD150136 


05 111(f) 


2569 


1.35 


43000 


16.0 


3.98 


3700 


2.02 


5 Per 


07.5 III(n)(f) 


4512 


0.40 


36000 


11.0 


0.32 


2600 


0.33 


t Ori 


09 III 


599, 2420 


0.50 


34000 


17.8 


1.10 


2000 


0.92 


P Cru 


B0.5 III 


2575 


0.15 


27500 


13.0 


0.05 


1600 


0.07 


Main Sequence 


















9Sgr 


04 V((f)) 


5398, 6285 


1.60 


43000 


16.0 


2.40 


2950 


1.53 


HD206267 


06.5 V((f)) 


1888, 189 


0.75 


40500 


12.7 


0.63 


3225 


0.46 


15 Men 


07V((f)) 


5401, 6248, 6243 


0.69 


40100 


9.7 


0.50 


2300 


0.67 


Ori C 


07 Vp 


3, 4 


0.55 


38000 


9.0 


0.20 


1650 


0.40 


COph 


09.5 Ve 


2571, 4367 


0.15 


34000 


8.0 


0.13 


1500 


0.30 


a Ori 


09.5 V 


3738 


0.50 


33000 


9.0 


0.08 


1250 


0.21 


T ScO 


BO V 


638, 2305 


0.17 


32000 


6.2 


0.03 


2400 


0.06 



Note. — The majority of stellar parameters are taken from Howarth & Prinja (1989) and Lamers & Leitherer (1993), along 
with input from Koch & Hrivnak (1981), Leitherer (1988), Howarth et al. (1993), CassineUi et al. (1994), and Waldron et al. 
(2004). Distances are taken from Savage et al. (1977), ShuU & Van Steenberg (1985), and Berghofer et al. (1996). 

Note. — Nwo is the scale factor of the stellar wind column density defined as Ny^o = M /i'K^ifj'mH'VooR* where iih and 
mn are respectively the mean molecular weight and mass of H atom. 
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Table 2. Lines Used in Analysis, Wind and ISM Cross Sections, and Relevant f /i 

Parameters 



ion 


Ao 


Tl 


aw 


criSM 




Nc 


J, 

<PC 




A 


MK 


1 n— 22 
10 


cm? 


J = 0, 1, 2 




photons 


S XVI 


4 727 4 733 


25.12 


0.205 


0.224 








S xv(r) 


5.039 


15.85 


0.240 


0.247 








S xv(i) 


5.063, 5,066 


12.59 


0.242 


0.249 








S xv(f) 


5.102 


15.85 


0.247 


0.255 


673.9, 738.2, 756.0 


1.9a;10" 


9.2x10^ 


Si XIV 


6.180, 6.186 


15.85 


0.401 


0.433 








Si xiii(r) 


6.648 


10.00 


0.424 


0.529 








Si xiii(i) 


6.685, 6.882 


10.00 


0.429 


0.536 








Si xiii(f) 


6.740 


10.00 


0.438 


0.496 


815.2, 865.2, 878.4 


4.0x10" 


2.4x10^ 


Mg XII 


8.419, 8.425 


10.00 


0.771 


0.912 








Mg xi(r) 


9.169 


6.31 


0.897 


1.148 








Mg xi(i) 


9.228, 9.231 


6.31 


0.912 


1.169 








Mg xi(f) 


9.314 


6.31 


0.933 


1.196 


997.7, 1034.3, 1043.3 


6.2x1012 


4.9x10^^ 


Ne X 


12.132, 12.138 


6.31 


1.825 


2.290 








Ne IX (r) 


13.447 


3.98 


2.138 


3.021 








Ne ix(i) 


13.550, 13.553 


3.98 


2.125 


3.082 








Ne ix(f) 


13.669 


3.98 


2.186 


3.174 


1247.8, 1273.2, 1277.7 


6.4x10" 


7.7x10^ 


Fe XVII 


15.014 


5.01 


2.193 


3.469 








Fe XVII 


16.780 


5.01 


2.917 


4.672 








VIII 


18.967, 18.972 


3.16 


3.999 


5.906 








Vll(r) 


21.602 


2.00 


1.735 


8.417 








Vll(i) 


21.801, 21.804 


2.00 


1.736 


8.631 








vii(f) 


22.098 


2.00 


1.796 


8.959 


1623.9, 1634.0, 1638.5 


3.4x10^° 


7.3x10^ 


N VII 


24.779, 24.785 


2.00 


2.387 


7.348 









Note. — Ao and Tl are taken from the APED line list. The A/_t are the EUV/UV wavelengths needed 

for the radiative excitation from the 2"^ 5*1 I^Pj levels (J = 0, 1, 2). Nc and (pc are the critical density 
and critical photo-excitation rate. If = then Nc represents the density required to reduce the f /i ratio 
by one- half, whereas, if the density is << Nc then (j)c represents the radiative excitation rate required to 
reduce the f /i ratio by one- half. 



Table 3. Observed MEG f /i Line Ratios and /ir-inferred Radii 



star 


O VII 




Ne IX 






Mg XI 












Si 


XIII 


S XV 




in 


Rfir 










/A 














Sli 






Rfir 


f lA 

J n 


Rfir 


Supergiants (I, II) 








































C Pup 


0.02 ± 0.02 


7.67 ± 3.29 


0. 


.35 ± 0.05 


10.56 ± 0.81 


0. 


33 it 0. 


,03 


3, 


,73 it 


0, 


,20 


1 


.16 ± 0. 


11 


2, 


,54 ± 0.21 


0.72 ± 0.24 


< 1.22 


Cyg OB2 No. 9 





























.61 ± 0. 


15 


1. 


,78 ± 0.27 






Cyg OB2 No. 8a 






0. 


,61 ± 0.24 


11.14 ± 3.26 


0. 


,86 ± 0. 


,22 


5. 


,61 ± 


0. 


,99 


1 


.02 it 0. 


13 


1. 


,91 ± 0.18 


0.75 ± 0.23 


< 1.18 


S Ori A 


0.10 ± 0.03 


13.32 ± 2.36 




< 0.001 


< 1.01 


0. 


81 ± 0, 


,19 


4, 


,43 ± 


0, 


,69 


1 


.44 ± 0. 


,42 


1, 


53 ± 0.43 






C Ori A 


0.12 ± 0.02 


13.51 ± I.IS 





27 ± 0.05 


6.15 ± 0.56 


1, 


10 ± 


14 


4, 


,97 ± 


0, 


,49 


2 


.62 ± 0. 


56 




> 1.62 


0.31 ± 0.26 


< 1.02 


c Ori 


0.10 ± 0.03 


10.14 ± 1.44 


0, 


51 ± 0.08 


7.60 ± 0.72 


1. 


21 ± 


24 


4, 


,47 ± 


0, 


73 


1 


.97 ± 0. 


52 


2 


96 ±1.96 






Giants (IV, III) 








































HD150136 






0, 


.12 ± 0.04 


5.86 ± 1.02 


0. 


,46 ± 0, 


,06 


4, 


,51 ± 


0, 


,37 


1 


.24 ± 0. 


17 


2, 


,81 ± 0.37 


2.50 ± 1.13 


> 1.76 


J Per 


< 0.01 


< 4 


0, 


.09 ± 0.03 


4.10 ± 0.74 


0. 


27 ± 0, 


.06 


2. 


.68 ± 


0, 


.28 


2, 


.61 ± 0. 


63 




> 3.44 






I Ori 


0.04 ± 0.02 


8.63 ± 2.14 


0, 


.62 ± 0.19 


11.14 ± 2.13 


0. 


,74 ± 0, 


.17 


4, 


.33 ± 


0, 


.67 


1 


.77 ± 0. 


,49 


2, 


,80 ± 1.32 






13 Cru 


< 0.03 


< 5 








0. 


28 ± 0, 


,17 


1, 


70 ± 


0, 


,55 


1 


.81 ± 1. 


13 




> 1.00 






Main Sequence 








































9 Sgr 


0.30 ± 0.12 


31.67 ± 7.03 


0, 


,02 ± 0.02 


2.49 ± 1.06 


0. 


63 ± 0, 


17 


.5 


46 ± 





92 


1 


.39 ± 0. 


37 


3, 


,31 ± 0.97 






HD206267 


< 0.03 


< 8.51 


0, 


,48 ± 0.19 


12.02 ± 2.94 


0. 


30 ± 


12 


3 


32 ± 





7fi 





.33 ± 0. 


18 




< 1.32 






15 Mon 


0.07 ± 0.05 


13.64 ± 5.32 




< 0.02 


< 2.15 


0. 


,29 ± 0, 


,24 


2, 


,96 ± 


1, 


,59 


1 


.79 ± 1. 


,44 




> 1.11 






Ori C 






0, 


.23 ± 0.08 


7.09 ± 1.31 


0. 


,12 ±0. 


.03 


1. 


.80 ± 


0. 


.21 


1 


.73 ± 0. 


16 


3. 


.26 ± 0.49 


1.38 ± 0.30 


1.53 ± 0.43 


C Oph 


0.00 


1.00 




< 0.02 


< 2.12 


0. 


,25 ± 0, 


,06 


2, 


,31 ± 


0, 


,28 


1 


.47 ± 0. 


,29 


1, 


,76 ± 0.38 


0.64 ± 0.41 


< 1.02 


(T Ori 


0.1)1) 


1.1)1) 





13 i l).()5 


4,51) ± 0,87 




< 0,1)4 




< 1, 


1)4 





,34 ± 1), 


24 




< 1.02 






r Sru 


O.Ul) 


l.Ol) 




1 O.Oi 




I), 


- 


05 




(ji = 





21) 




,4!) ± 1), 






2.:!,"'. 


i.7,-, ± (),,-,7 


1 i.OO 



I 



Table 4. Observed HEG f /i Line Ratios and /ir-inf erred Radii 



Star 


Ne IX 


Mg XI 


Si XIII 


S XV 






f/i 




f/i 




f/i 




f/i 


Rfir 


Supergiants (I, II J 




















0.13 ±0.05 


6.15 ±1.16 


0.35 ± 0.06 


3.79 ± 0.39 


0.98 ±0.16 


2.20 ±0.28 


1 1 Q _i_ n Rn 


Z.lo 


Cyg (Jr52 JNo. 9 






< 1.30 


< 8.81 


< 0.99 


< 2.06 






/"I *^T!>ri AT— 0_ 

Cyg (Jh!2 INo. »a 


< 0.09 


< 4.33 


0.96 ± 0.34 


6.12 ± 1.58 


1.15 ±0.20 


2.11 ±0.33 


U.oo ± U.z9 


< 1.20 


Uri A 






0.86 ± 0.47 


4.61 ± 1.77 


2.54 ±1.30 


> 1.23 






C, UTl A 


0.30 ±0.10 


6.44 ± 1.20 


0.62 ±0.18 


3.37 ± 0.60 


1.51 ±0.52 


< 1.56 


< 0.16 


< 1.02 


^ Oi-; 
e Uri 


0.30 ±0.12 


5.56 ±1.21 


2.12 ±0.64 


11.27 ±5.94 


3.71 ±2.35 


> 1.00 






Giants (IV, 111 J 


















HD150136 


0.20 ±0.13 


7.47 ± 2.82 


0.35 ± 0.08 


3.87 ±0.51 


1.16 ±0.27 


2.67 ±0.55 






^ Per 






0.33 ±0.10 


2.89 ± 0.50 


0.31 ±0.11 


< 1.02 






L Ori 


0.55 ± 0.43 


9.86 ± 5.38 


0.69 ± 0.20 


4.14 ±0.78 


< 0.02 


< 1.02 






/3 Cru 


















Main Sequence 


















9Sgr 


< 0.11 


< 5.85 


0.54 ±0.22 


4.98 ± 1.24 










HD206267 


















15 Mon 


< 0.20 


< 7.20 


0.78 ± 0.42 


5.84 ± 2.20 










el Ori C 


< 0.01 


< 1.20 


0.41 ± 0.10 


3.51 ± 0.52 


1.10 ±0.14 


1.96 ±0.21 


0.68 ±0.17 


< 1.02 


COph 


0.18 ±0.11 


5.38 ± 1.87 


0.27 ±0.10 


2.38 ± 0.49 


1.37 ±0.49 


1.76 ±0.61 






CT Ori 


0.16 ±0.15 


4.38 ±2.87 














T ScO 


0.05 ± 0.03 


2.53 ± 0.98 


0.29 ±0.10 


2.31 ±0.41 


1.77 ±0.34 


1.69 ±0.49 







Table 5. Observed MEG H/He Line Ratios and Derived Thhc 



Star 




Oxyg 


en 






Neon 






Magnesium 






Sil 


con 


Sulfur 






H/He 










H /He 








H/He 






H/He 






H/He 




Tl range 






2.00 - 3.16 






3.98 - 6.31 






6.31 - 10.00 








10.00 - 16.86 




15.85 - 25.12 


Supergiants (I, II) 






































C Pup 


1 


57 ± 0.20 


2 


75 ± 


10 





84 ± 0.05 


3.81 ± 


06 





40 ± 0.03 


6.07 ± 0.12 





16 ± 


02 


7.68 ± 0.22 






Cyg OB2 No. 9 




















3 


26 ± 0.92 


12.48 lb 1.53 


1 


28 ± 


30 


13.93 ± 1.25 






Cyg OB2 No. 8a 




















1 


23 ± 0.15 


8.50 ± 0.38 


1 


08 it 


09 


13.19 it 0.42 


1.60 ± 0.38 


23.05 ± 2.42 


6 Ori A 


1 


34 ± 0.13 


2 


62 ± 


07 





80 ± 0.09 


3.76 ± 


11 





31 ± 0.05 


5.73 it 0.24 





57 ± 


14 


10.64 ± 0.80 






C Ori A 


1 


14 ± 0.07 


2 


51 ± 


04 





48 ± 0.04 


3.28 ± 


07 





26 ± 0.04 


5.50 ± 0.18 





21 it 


05 


8.08 it 0.53 






e Ori 


1 


18 ± 0.10 


2 


53 ± 


06 





57 ± 0.05 


3.42 ± 


09 





17 ± 0.03 


5.01 ± 0.19 





03 ± 


02 


5.07 ± 0.97 






mean T/j^g 






2 


60 ± 


07 






3.57 ± 


08 






7.22 ± 0.66 








9.76 ± 0.78 




23.05 ± 2.42 


Giants (IV, III) 






































HD150136 












1 


68 ± 0.18 


4.68 ± 


15 





44 ± 0.04 


6.20 ± 0.14 





50 ±0 


06 


10.27 ± 0.36 


0.76 ± 0.25 


17.71 ± 2.10 


5 Per 


2 


74 ± 0.43 


3 


26 ± 


17 





92 ± 0.09 


3.91 ± 


11 





16 ± 0.03 


4.87 ± 0.17 





11 ± 


04 


6.99 ± 0.64 






1, Ori 


1 


35 ± 0.13 


2 


63 ± 


OS 





18 ± 0.03 


2.58 ± 


08 





32 ± 0.08 


5.73 ± 0.36 





64 ±0 


20 


11.00 ± 1.09 






13 Cru 


1 


19 ± 0.18 


2 


54 ± 


11 





44 ± 0.09 


3.19 ± 


17 




















mean ThHb 






2 


81 ± 


13 






3.69 ± 


13 






5.60 ± 0.24 








9.42 ± 0.63 




17.71 ± 2.10 


Main Sequence 






































9 Sgr 


2 


63 ± O.BB 


3 


22 ± 


22 





45 ± 0.06 


3.22 ± 


10 





30 ± 0.05 


5.68 ±0.24 





17 ± 


06 


7.72 ±0.60 






HD 206267 












1 


36 ± 0.27 


4.38 ± 


26 





24 ± 0.08 


5.36 ± 0.43 





13 ± 


08 


7.06 ± 1.06 






15 Mon 





74 ± 0.11 


2 


24 ± 


08 





12 ± 0.04 


2.41 ± 


14 





08 ± 0.05 


4.25 ± 0.54 














Ori C 












2 


46 ± 0.25 


5.32 ± 


21 


2 


36 ± 0.18 


10.97 ± 0.35 


1 


76 ± 


09 


15.80 ± 0.31 


1.39 ± 0.16 


22.49 ± 1.06 


C Oph 


1 


73 ± 0.28 


2 


83 ± 


14 


1 


39 ± 0.14 


4.42 ± 


13 





51 ± 0.08 


6.49 ± 0.28 





38 ± 


08 


9.48 ± 0.50 






(T Ori 


1 


76 ± 0.32 


2 


84 ± 


16 





69 ± 0.12 


3.60 ± 


17 





04 ± 0.03 


2.13 ± 2.12 














T ScO 


3 


81 ± 0.49 


3 


66 ± 


17 


1 


73 ± 0.12 


4.72 ± 


10 





84 ± 0.07 


7.50 ± 0.19 





61 it 


06 


10.92 ± 0.35 


0.86 it 0.26 


18.57 it 2.07 


mean Thh^ 






2 


96 ± 


16 






3.79 ± 


16 






5.23 ± 0.86 








8.80 ± 0.66 




18.57 ± 2.07 


Note. — H/He rat 


OS are ISM corrected and T/f 


H 


g are the derived X-ray temperatures in 


MK. The results for 


Or 


C are not included 


in the mean T/f 


He ■ 



Table 6. Observed HEG H/He Line Ratios and Derived Thhc 



Star Neon Magnesium Silicon Sulfur 





H/He 


He 


H/He 


He 


H/He 


^H He 


H/He 


'^H He 


Tl range 




3.98 - 6.31 




6.31 - 10.00 




10.00 - 15.85 




15.85 - 25.12 


Supergiants (I, II) 


















C Pup 


0.80 ± 0.09 


3.76 ±0.12 


0.28 ± 0.04 


5.56 ±0.18 


0.10 ±0.03 


6.86 ± 0.41 






Cyg OB2 No. 9 


















Cyg OB2 No. 8a 






1.30 ±0.24 


8.67 ±0.60 


1.08 ±0.16 


13.16 ±0.69 


0.82 ± 0.27 


18.25 ±2.21 


<5 Ori A 


0.18 ± 0.06 


2.60 ±0.19 


1.00 ±0.33 


7.90 ±0.87 


0.24 ±0.10 


8.34 ± 0.92 






C Ori A 


0.39 ± 0.06 


3.10 ±0.12 


0.28 ±0.06 


5.55 ± 0.29 


0.30 ±0.10 


8.81 ±0.81 






e Ori 


0.48 ± 0.09 


3.28 ±0.16 


0.26 ±0.07 


5.47 ± 0.37 


0.10 ± 0.05 


6.81 ± 0.80 






mean Tjju^ 




3.18 ±0.15 




6.63 ± 0.52 




8.80 ± 0.75 




18.25 ±2.21 


Giants (IV, III) 


















HD150136 


0.54 ±0.16 


3.36 ±0.27 


0.46 ± 0.08 


6.29 ± 0.28 


0.50 ±0.09 


10.28 ±0.59 






e, Per 


0.34 ± 0.09 


3.00 ± 0.19 


0.22 ± 0.07 


5.27 ±0.39 


0.23 ±0.10 


8.23 ± 0.95 






(, Ori 


0.86 ± 0.23 


3.81 ±0.29 














P Cru 


















mean Thhr 




3.39 ± 0.25 




6.00 ± 0.34 




9.26 ± 0.79 






Main Sequence 


















9 Sgr 


1.06 ± 0.27 


4.06 ± 0.31 


0.22 ± 0.08 


5.23 ± 0.42 


0.37 ±0.13 


9.38 ± 0.94 






HD206267 


















15 Mon 


0.27 ±0.11 


2.80 ± 0.29 














6li Ori C 






3.07 ± 0.38 


12.25 ±0.64 


1.68 ±0.13 


15.52 ±0.46 


1.00 ±0.15 


19.74 ± 1.09 


COph 


1.14 ±0.22 


4.15 ±0.25 


0.69 ±0.15 


7.04 ± 0.47 


0.29 ±0.10 


8.74 ± 0.81 






T ScO 


1.79 ± 0.26 


4.76 ±0.21 


0.84 ±0.12 


7.49 ± 0.32 


0.68 ±0.10 


11.28 ±0.55 


0.15 ±0.09 


10.89 ± 1.71 


mean TnHe 




3.94 ±0.27 




6.59 ± 0.41 




9.80 ± 0.78 




10.89 ±1.71 


Note. — H/He ratios are ISM corrected and Thhb are the derived X-ray temperatures in 


MK. The results for 6^ Ori C 


are not included 



in the mean Tuu^. 



Table 7. Observed MEG and HEG G-Ratios 



Star 




Oxygen 
MEG HEG 




MEG 


Neon 


HEG 






Magnesium 
MEG HEG 






MEG 


Sil 


con 


HEG 






MEG 


Sulfur 

HEG 


Siipc.rqiants (I, II) 
















































Q Pup 


1 


09 ±0.16 





48 ± 


04 





64 ± 


09 





82 ± 0.05 





77 ± 


08 





96 ± 


06 


1 


14 ± 


14 


2 


45 ± 


67 


1.21 ± 0.47 


Cyg OB2 No. 9 





















50 ± 0.21 








1 


90 ± 


40 
















Cyg OB2 No. 8a 





















40 ± 0.06 





68 ± 


15 





90 ± 


08 


1 


31 ± 


18 


1 


87 ± 


46 


1.99 ± 0.56 


5 Ori A 





79 ± 0.09 





60 ± 


08 





44 ± 


16 





62 ± 0.09 





84 ± 


30 





91 ± 


18 





67 ± 


19 


2 


11 ± 1 


51 


0.75 ± 0.50 


C Ori A 





94 ± 0.07 





64 ± 


06 





64 ±0 


09 


1 


03 ± 0.09 





69 ± 


12 


1 


21 ± 


17 





92 ±0 


22 


1 


43 ± 


80 


< 0.48 


e Ori 





98 ±0.10 





86 ±0 


08 





62 ± 


12 





75 ± 0.09 





97 ± 


19 


1 


02 ± 


18 





43 ± 


13 










Giants (IV, III) 
















































HD160136 









87 ± 


11 


1 


01 ± 


32 





64 ± 0.05 





95 ± 


13 





76 ± 


07 





88 ± 


14 


1 


20 ± 


36 


0.70 ± 0.32 


e Per 





53 ±0.13 





66 ± 


08 





77 ± 


21 





79 ± 0.09 


1 


19 ± 


23 





91 ± 


13 


1 


33 ± 


34 


1 


02 ± 


66 




L Ori 





91 ±0.11 





18 ± 


03 





31 ± 


13 





90 ± 0.14 


3 


52 ± 1 


05 


1 


85 ± 


42 


1 


06 ± 


40 










(3 Cru 


1 


44 ± 0.24 





96 ± 


18 








1 


03 ± 0.37 








1 


34 ± 


61 
















Main Sequence 
















































9 Sgr 


1 


46 ± 0.35 





45 ±0 


06 





71 ±0 


23 





32 ± 0.05 





77 ± 


19 





52 ± 


09 





67 ±0 


19 










HD206267 





18 ±0.14 





69 ± 


16 











70 ± 0.15 











66 ± 


20 
















15 Mon 





57 ±0.11 





94 ± 


18 





53 ± 


22 





56 ± 0.23 


2 


95 ± 1 


59 





68 ± 


33 





62 ± 


40 










el Ori C 


1 


40 ± 0.49 





51 ± 


08 





89 ± 


34 





69 ± 0.07 





95 ± 


15 





72 ± 


04 





92 ± 


08 





73 ± 


10 


0.83 ± 0.13 


C Opli 





92 ±0.19 





67 ± 


08 





60 ± 


16 





78 ± 0.09 


1 


00 ± 


22 





98 ± 


13 





84 ± 


20 


1 


99 ± 1 


06 




CT Ori 


1 


09 ± 0.24 


1 


12 ± 


18 





58 ± 


22 





89 ± 0.23 


1 


06 ± 


53 


1 


26 ± 


59 





74 ± 


41 










r Sco 


1 


59 ± 0.25 





89 ± 


08 





78 ± 


15 





80 ± 0.07 





58 ± 


10 


1 


04 ± 


09 


1 


14 ± 


15 


1 


80 ± 


47 


0.53 ± 0.22 



00 

Note. — G-Ratios are ISM corrected. ^ 
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Fig. 1.— The MEG and HEG HWHM, Vp, and EMx histograms for all OB stars illustrating 
the percentage of lines within a given parameter bin range. 
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Fig. 2. — The MEG and HEG HWHM /voo histograms illustrating the luminosity class 
dependence of the percentage of lines within a given bin range. 
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0.6 




Vp/V» 



Vp/V» 




Fig. 3. — The MEG and HEG Vp/v^o histograms illustrating the luminosity class dependence 
of the percentage of lines within a given bin range. 
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Fig. 4. — The MEG and HEG log EMx histograms illustrating the luminosity class depen- 
dence of the percentage of lines within a given bin range. The EMx have units of cm~^. 
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X-ray Temperature (MK) 

Fig. 5. — The H/He line ratio dependence on Tx for O, Ne, Mg, Si, and S as determined 
from the APED data. For the He-hke hnes we only use r-line (see text). For each H-like and 
He-like line region all lines within the instrumental wavelength resolution are included in 
determining the corresponding total emission at a given temperature. Since the H-like Ne x 
line (~ 12.13 A) is a blend with a relatively strong Fe xvil (~ 12.12 A) and is unresolvable 
with the HETGS, the dashed line shows the expected temperature dependence of the Ne 
H-He ratio if the lines were resolved, whereas, the solid line represents the realistic observed 
ratio (i.e., the ratio of the Fe xvil and Ne x line emission to the Ne ix r-line emission). For 
Tx > 8 MK, the Ne ix r-line is enhanced by other Fe lines producing the drop in the Ne ix 
H/He ratio. 
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Fig. 6. — MEG and HEG scatter plots showing the relation between Rfij. and Rt=i for all 
luminosity classes. The error bars for Rr=i are determined for a X-ray continuum optical 
depth of 1 ± 0.5 and the adopted value of Rt=i is the average of these limits. 
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Fig. 7. — Scatter plots showing the dependence of the MEG derived He-hke ion HWHM 
(top panel) and Vp (bottom panel) on the ambient wind velocity, Vq, for all OB luminosity 
classes. All quantities are determined for each at their Rfij. value and normalized to their 
star's Voo- The dashed-lines represent special cases when HWHM = Vq and Vp = —Vq 
(negative means the expected maximum blue-shift velocity) . The HWHM range from to 
0.6 with no clear indication of any dependence on Vq. The HWHM < 0.2 v^o pose a 
problem since they are > the local ambient velocities at these low radial positions. All Vp 
are concentrated about zero velocity with no indication of any large blue-shifted lines. 
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Fig. 8. — MEG and HEG scatter plots showing the dependence of Tx (MK) (determined 
from H/He line ratios) on the He-hke /ir-inferred radii, Rfir (determined from He-hke f/i 
hne ratios), for all OB luminosity classes. 
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Fig. 9. — Scatter plot of U^ei (pre-shock velocity relative to the shock front) versus Vq (the 
ambient wind velocity) for all OB luminosity classes. The values of Vrd are determined from 
the MEG derived values of Tuue using eq. [5l Both Urei and Vq for each star are evaluated 
at their Rji^ value and normalized to their star's foo- The dashed-lines represent two special 
cases where U^ei = and [/^e/ = 0.5Vo. A detailed discussion is given in Section 6. 
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A. ERRATUM: "An Extensive Collection of Stellar Wind X-ray Source 

Region Emission Line Parameters, Temperatures, Velocities, and Their 
Radial Distributions as Obtained from Chandra Observations of 17 OB 

Stars" (ApJ, 608, 456 [2007] ) 

The major objective of the paper was to provide a detailed tabulation of the observed 
HETGS X-ray emission hne flux ratios. We presented the MEG and HEG He-like f /i line 
ratios, the H-hke to He-hke {H/ He) hne ratios, and the He-hke G-ratios. The stellar wind 
spatial locations of the X-ray sources were derived from the f /i ratios and their associated 
X-ray temperatures were obtained from the H/He ratios (Thhc)- This information was used 
to verify the correlations between Rfi^ and Rr=i (Figure 6) and T^He and Rfir (Figure 8). 
However, we have realized that some of our tabulated uncertainties for these line ratios were 
underestimated, primarily for those lines with low S/N data. Hence, the primary purpose 
of this erratum is to provide a tabulation of the corrected line ratios and their uncertainties. 

The details of our line fitting procedure are discussed in Section 3.3. All uncertainties 
were determined using standard statistics (e.g., Bevington 1969). First, we would like 
to clarify a statement in Section 3.3 (second paragraph), which states that all parameter 
uncertainties were determined from 90% confidence regions, but in actuality all uncertainties 
were estabhshed using 68% confidence regions. With regards to the main point of this 
erratum, we found that our algorithm for determining the covariance matrix which is 
used to determine the uncertainties of the fitting parameters had an indexing error in the 
coding logic which produced errors in some of the off-diagonal terms. From our detailed 
examination of the code, we found that certain cases were especially vulnerable to this 
coding error, in particular, those cases where the normalization ranges were large (i.e., 
low S/N data). This code correction has also produced changes in some of the line ratios 
and their derived quantities (e.g., Rfir and Thhc)- 

The algorithm has been corrected and the aff'ccted Tables (Tables 3, 4, 5, 6, and 7) have 
been updated and re-produced in this erratum. We also corrected a few entries that were 
originally tabulated incorrectly, and some data were removed as they did not satisfy our S/N 
criterion, i.e., the HEG S XV f/i data for ( Ori, the MEG S XV f/i data for ( Oph, and the 
HEG Mg XI and Si xill f/i data for Cyg 0B2 No. 9. As discussed in Section 3.1, we stated 
that if a reasonable fiux had been estabhshed, these results would be used only for estimating 
hne ratios that provide interesting limits. However, the meaning of a "reasonable" fiux limit 
was unclear; the criterion used is that the observed total net counts from all three He-like fir 
lines must have a S/N > 3. We also need to clarify the significance of blank entries in our 
Tables. The blanks just indicate that the given line ratio has either an unphysical result that 
produces an anomalously large uncertainty, or did not satisfy our He-like S/N >3 criterion. 
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An unphysical result example is when the fitting procedure predicts a He-like i-line flux that 
is too small. This occurs primarily in low S/N high energy He-like fir lines where the effects 
of line overlap can lead to a poor determination of the i-line. 

In addition, we would like to clarify why the relative uncertainties in Rfir are typically 
smaller than the corresponding f /i relative uncertainties. As shown in eq. (2) of the paper, 

for (f)/(j)c > 1 si'iid Tie << Nq (valid throughout the wind except when extremely close to 
the star), the f /i ratio is inversely proportional to the dilution factor, W(r), which changes 
rapidly for small changes in radius. Hence, dramatic changes in the f /i ratio can occur for 
only small changes in radius, which explains the differences seen in the relative uncertainties. 
Plus, another key point that was not mentioned in the original paper is that there is a lower 
hmit on the f /i ratio for the case where Ue « Nq determined by setting W = 0.5. This 
implies that any observed f/i ratio below this limit indicates that density effects may be 
important. We are currently investigating this possibility. 

The format of the data presented in these new tables have been changed slightly: 1) 
all best-fit line ratios {f/i and H/He) and their uncertainties are tabulated, regardless of 
the size of the uncertainty; 2) as before, all derived Rfir and Thhc represent an average 
of their respective ranges predicted by the f/i and H/He line ratios, and their associated 

uncertainties arc equal to half the difference in these ranges; 3) for those cases where the f/i 
uncertainty is > the best-fit f/i ratio, the derived Rfir average is determined using a lower 
limit of Rfif = i?*, and an upper limit on Rfir determined by the f/i + uncertainty; 5) for 
those cases where the minimum f/i {f/i - uncertainty) ratio predicts a finite Rfir, but the 
upper limit in the f/i range is at or greater than its asymptotic value (i.e., the low-density 
and zero UV flux limit), the upper limit on Rfir is undetermined (i.e., Rfir oo), and these 
Rfir values are tabulated as lower limits, and; 5) for those cases where the H/He uncertainty 
is > the H/He ratio, upper limits on Thhc a-re presented. 

The two key figures of the paper (Figures 6 and 8) have also been re-produced. In 
these figures only data with finite limits on Rfir and Thhb are plotted, i.e., data with 
just lower limits on Rfir and upper limits on Thhc are not shown. For clarity, we also 
chose not to display any Rfir data where the uncertainty is > lOi?*. The impact of these 
changes in the other figures that depend on the new derived Rfir and uncertainties (Figure 
7 and 9) are found to show minimal differences from the original results. However, we 
did find an erroneous high temperature data point in Figure 9 at Vo{Rfir)/voo ~ and 
Urei{TH He) /voo ~ 0.82 and it should be ignored. This same data point at low Rfir was also 
in the original Figure 8 for the giants in both the MEG and HEG plots, and it has been 
removed from the re-produced Figure 8. The source of this data point was traced to the 
star 7 Vel, originally considered in our analysis, but was dropped from our study due to its 
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highly unusual X-ray spectra which were deemed inappropriate for this study of "normal" 
OB stars. We have confirmed that no other data points from 7 Vel were present in any of 
the original plots. 

We would also like to add a comment concerning the importance of obtaining high 
S/N HEG data as illustrated by comparing the S XV MEG and HEG determined f/i ratios 
for 0^ Ori C. This is a clear example of how line overlap, either caused by the physical line 
width or the energy resolution capabilities of the instrument, can lead to larger uncertainties. 
Although for this case both the MEG and HEG do have comparable S/N data, the MEG data 
had significant line overlap produced by the best-fit line width and lower energy resolution 
of the MEG, whereas, in the HEG data the S xv fir were resolved leading to a significant 
reduction in the uncertainty. 

In general, comparisons of these new derived quantities with the original tabulated data 
show that the largest differences are seen in the uncertainties, primarily the uncertainty 
results for O Vll and S XV and any other low S/N fines. In addition, there are minor 
changes evident in some of the f/i ratios, the H/He ratios, and the G-ratios. In addition, 
the newly tabulated data and the reproduced Figures 6 and 8 indicate that the fundamental 
results discussed in the paper have not changed, and these corrections have not altered any 
of the conclusions discussed in the paper. 

We wish to thank Maurice Leutenegger for his communication concerning the specific 
details of our line fitting approach. This inquiry motivated us to re-examine all aspects of 
our line fitting algorithms whereupon we found the above mentioned coding mistake in the 
line fiux uncertainties. 

Bevington, P. R. 1969, "Data Reduction and Error Analysis for the Physical Sciences" , 
McGraw-Hill, Inc. 



Table 3. Observed MEG f /i Line Ratios and /ir-inferred Radii 



star 




O VII 




Ne IX 




Mg XI 


Si XIII 


S XV 




f/i 


Rfir 


f/i 


Rfir 


f/i 


Rfir 


f/i Rfir 


f/i Rfir 



Supergiants (I, II) 



C Pup 





03 


± 





11 


11.21 ± 10.21 





38 


± 





13 


10.85 ± 2.15 





34 


± 





05 


3 


76 


± 


0.30 


1 


17 


± 





17 


2.56 ± 0.33 


0.60 ± 0.62 


1.24 ± 0.24 


Cyg OB2 No. 9 















































60 


± 





20 


1.49 ± 0.28 






Cyg OB2 No. 8a 

















51 


± 


1 


29 


16.32 ± 14.32 





79 


± 





65 


6 


22 


± 


3.20 


1 


02 


± 





22 


1.91 ± 0.31 


0.74 ± 0.34 


1.07 ± 0.07 


S Ori A 





11 


± 





08 


12.98 ± 5.45 





00 


± 





02 


1.43 ± 0.43 





84 


± 





27 


4 


49 


± 


0.98 


1 


45 


± 





54 


1.63 ± 0.59 






C Ori A 





13 


± 





03 


13.85 ± 1.60 





32 


± 





10 


6.76 ± 1.16 


1 


10 


± 





18 


4 


97 


± 


0.64 


3 


48 


± 


1 


34 


> 1.80 


0.77 ± 1.21 


1.33 ± 0.33 


c Ori 





10 


± 





05 


10.23 ± 2.56 





51 


± 





11 


7.54 ± 0.99 


1 


24 


± 





35 


4 


61 


± 


1.09 


1 


98 


± 





68 


> 1.0 






Giants (IV, III) 




























































HD150136 





00 


± 





10 


9.80 ± 8.80 





12 


± 





10 


5.30 ± 2.96 





45 


± 





16 


4 


41 


± 


0.90 


1 


26 


± 





40 


2.97 ± 0.91 


1.46 ± 2.25 


> 1.0 


J Per 





00 


± 





09 


7.62 ± 6.62 





09 


± 





08 


3.69 ± 2.19 





26 


± 





10 


2 


52 


± 


0.62 


2 


36 


± 


1 


06 


> 1.90 






I Ori 





03 


± 





08 


7.64 ± 6.64 





58 


± 


1 


06 


12.29 ± 11.29 





72 


± 





23 


4 


24 


± 


0.90 


1 


79 


± 





90 


> 1.16 






/3 Cru 





01 


± 





04 


4.03 ± 3.03 





00 


± 





06 


1.67 ± 0.67 





32 


± 





26 


1 


79 


± 


0.79 


1 


43 


± 


1 


27 


> 1.0 






Main Sequence 




























































9 Sgr 





22 


± 





44 


25.75 ± 24.75 





05 


± 





11 


4.06 ± 3.06 





84 


± 


1 


96 


19 


14 


± 


18.14 


1 


43 


± 


1 


31 


> 1.0 






HD206267 





00 


± 





56 


22.12 ± 21.12 





56 


± 





45 


12.61 ± 7.08 





30 


± 





21 


3 


17 


± 


1.38 





34 


± 





32 


1.26 ± 0.26 






15 Mon 





05 


± 





26 


16.00 ± 15.00 





01 


± 





11 


3.35 it 2.35 





22 


± 





62 


3 


56 


± 


2.56 


1 


16 


± 


1 


85 


> 1.0 






Ori C 





00 


± 





12 


8.96 ± 7.96 





22 


± 





21 


5.92 ± 4.22 





12 


± 





09 


1 


72 


± 


0.72 


1 


73 


± 





27 


3.40 ± 0.84 


1.32 ± 0.66 


2.94 ± 1.94 


C Oph 





00 


± 





16 


9.25 ± 8.25 





01 


± 





05 


2.12 ± 1.12 





26 


± 





08 


2 


36 


± 


0.38 


1 


46 


± 





34 


1.76 ± 0.44 






(T Ori 


1) 


00 


± 





1)3 


3.92 ± 2.92 


1) 


13 


± 


1) 


06 


4.37 ± 1.19 





1)2 


± 


1) 


10 


1 


2!) 


± 


0.29 





21 


± 







<- 1.02 






r Sru 


11 


00 







IJi 




11 






11 


0:i 


l.Ui ^ 





:i6 




I) 


07 




">•> 




0.27 











40 


^ 1.7(j 


i.74 ± 1).,^4 


^ 1.0 



I 



Table 4. Observed HEG f /i Line Ratios and /ir-inf erred Radii 



Star 


Ne IX 


Mg XI 


Si XIII 


S XV 




f/i 






R( 




Rt 


f/i 


Rt 


Supergiants (I, II) 




















0.13 ±0.10 


5.74 ± 2.55 


0.36 ± 0.09 


3.84 ± 0.53 


0.97 ± 0.20 


2.19 ±0.35 


0.79 ± 0.65 


1.42 ± 0.42 


/~1 /^"Dri "NT— r\ 

Oyg UiiZ J\o. 9 


















Oyg (Jhi2 INo. 8a 






0.91 ±0.61 


5.95 ±2.97 


1.15 ±0.33 


2.15 ±0.53 


0.87 ±0.37 


1.13 ±0.13 


S Ori A 


0.00 ±0.22 


3.52 ± 2.52 


0.79 ±0.93 


4.64 ± 3.64 


2.52 ±1.91 


> 1.0 






C, Uri A 


0.31 ±0.17 


6.40 ±2.12 


0.62 ±0.24 


3.32 ±0.82 


1.82 ±1.11 


> 1.0 






e Uri 


0.30 ± 0.25 


5.09 ± 2.80 


2.08 ± 0.84 


15.73 ±11.24 










Lriants (IV, 111 J 


















HD150136 


0.27 ± 0.60 


9.77 ± 8.77 


0.34 ±0.12 


3.78 ± 0.77 


1.15 ±0.58 


2.85 ±1.27 






^ Per 


0.00 ± 0.09 


2.58 ± 1.58 


0.35 ±0.16 


2.92 ±0.78 


0.30 ±0.26 


1.06 ± 0.06 






1, Ori 


0.57 ±1.94 


22.84 ± 21.84 


0.68 ± 0.33 


4.10 ±1.30 


2.50 ±2.73 


> 1.0 






/3 Cru 


0.10 ±0.53 


4.70 ± 3.70 


0.55 ± 0.54 


2.47 ± 1.47 










Main Sequence 


















9Sgr 


0.05 ± 0.50 


7.40 ± 6.40 


0.54 ± 0.41 


4.76 ± 2.46 










HD206267 






1.02 ± 1.00 


7.56 ± 6.52 


1.97 ± 1.71 


> 1.02 






15 Mon 


0.00 ± 0.39 


5.83 ±4.83 


0.76 ± 0.94 


6.12 ±5.12 










01 Ori C 


0.00 ±0.14 


3.31 ± 2.31 


0.38 ±0.20 


3.24 ±1.01 


1.11 ±0.19 


1.97 ±0.28 


0.70 ± 0.25 


1.02 ± 0.02 


COph 


0.22 ± 0.33 


5.75 ±4.75 


0.27 ±0.14 


2.37 ±0.68 


1.37 ±0.57 


1.83 ±0.73 






a Ori 


0.17 ±0.43 


5.80 ±4.80 


0.00 ±0.12 


5.80 ±4.80 










T ScO 


0.03 ± 0.09 


2.59 ± 1.59 


0.29 ±0.13 


2.26 ±0.56 


1.77 ±0.37 


1.72 ±0.54 







Table 5. Observed MEG H/He Line Ratios and Derived Thhc 



Star 




Oxyg 

H/He 


en 




Neon 

H/He 








Magnesium 
H/He Thhc 




H/He 


Sili 


icon 

THHe 




Sulfur 
H/He THHe 


Tl range 






2.00 - 3.16 






3 


.98 - 6.31 








6.31 - 10.00 








10.00 - 16.86 






15.85 - 25.12 


Supergiants (I, II) 










































C Pup 


1 


.50 ± 0.27 


2, 


,68 ± 0.13 





,96 ± 0.10 


3. 


,92 ± 0, 


,12 


0, 


,37 ± 


0.05 


5.86 ± 0.20 


0. 


16 ± 0, 


,03 


7.66 ± 0.42 


0. 


,09 ± 0.15 


< 13 


Cyg OB2 No. 9 




















0. 


,82 ± 


0.53 


7.12 ± 1.61 


1. 


01 ± 0, 


,33 


12.63 ± 1.49 








Cyg OB2 No. 8a 























91 ± 


0.22 


7.59 it 0.63 


0. 


91 it 0. 


12 


12.34 it 0.57 


2. 


25 it 1.44 


26.28 it 8. 12 


6 Ori A 


1 


,21 ± (J. 13 


2 


,53 ± 0,07 





.90 ± 0.13 


3 


84 ± 


17 


t) 


29 ± 


0.07 


5.53 ± 0.31 


0. 


58 ± 0. 


18 


10.57 ± 1.01 








C Ori A 





.99 ± 0.06 


2 


,38 ± 0.05 





.66 ± 0.06 


3, 


,52 ± 0, 


,09 


0. 


,25 ± 


0.05 


5.32 ± 0.23 


0. 


19 ± 0. 


,07 


7.79 it 0.70 


1, 


,30 ± 4.13 


< 41 


e Ori 


1 


.02 ± 0.09 


2 


.41 ± 0.07 





.72 ± 0.07 


3, 


,60 ± 0, 


,10 


0, 


,18 ± 


0.04 


4.96 ± 0.26 


0. 


03 ± 0, 


,10 


< 7.2 


0, 


,13 ± 1.07 


< 22 


mean T/j^g 






2 


.50 ± 0.11 






3, 


,72 ± 0, 


,16 








6.06 ± 0.83 








10.18 lb 1.71 






26.28 ± 8.12 


Giants (IV, III) 










































HD150136 


1, 


.94 ± 1.06 


2, 


.83 ± 0.51 


1, 


.27 ± 0.24 


4, 


,26 ± 0, 


,22 


0, 


,40 ± 


0.08 


6.01 ± 0.34 


0, 


,49 ± 0, 


,09 


10.11 ± 0.52 


1, 


,18 ± 1.41 


< 30 


5 Per 


1 


.49 ± 0.36 


2 


.67 ± 0.18 





.74 ± 0.11 


3, 


,63 ± 0, 


,16 


0, 


,15 ± 


0.05 


4.77 ± 0.34 


0. 


,14 ± 0, 


,09 


7.08 ± 1.11 








1, Ori 


1 


.27 ± 0.16 


2 


.56 ± 0.08 





.36 ± 0.13 


2 


,97 ± 


,29 


0, 


,28 ± 


0.10 


5.44 ± 0.48 


0, 


,67 ± 0, 


,32 


10.92 ± 1.82 








13 Cru 


1 


.12 ± 0.18 


2 


.47 ± 0.11 


0, 


.41 ± 0.10 


3, 


,09 ± 0, 


,20 


0, 


,21 ± 


0.26 


< 6.4 
















mean ThHb 






2, 


.63 ± 0.28 






3, 


,48 ± 0, 


,42 








5.41 ± 0.41 








9.37 ± 1.29 








Main Sequence 










































9 Sgr 


0, 


.84 ± 0.33 


2, 


.28 ± 0.24 


0, 


.43 ± 0.09 


3, 


,13 ± 0, 


,18 


0, 


,27 ± 


0.13 


5.38 ±0.66 


0, 


14 ±0, 


,07 


7.21 ±0.83 








HD 206267 


0, 


.47 ± 0.23 


1 


.94 ± 0.23 


0, 


.32 ± 0.09 


2, 


,91 ± 0, 


,22 


0, 


,14 ± 


0.08 


4.62 ± 0.67 


0, 


15 ± 0, 


,20 


< 9.2 








15 Mon 


0, 


.71 ± 0.15 


2 


.18 ± 0.11 


0, 


.18 ± 0.07 


2 


52 ± 0, 


21 





12 ± 


0.13 


< 5.4 


0, 


14 it 2, 


,37 


< 18 








fll Ori C 


3, 


.09 ± 1.84 


3, 


.28 ± 0.73 


2 


.53 ± 0.36 


.5 


33 ± 


26 


2 


,06 ± 


0.22 


10.31 ± 0.43 


1, 


68 ± 0, 


11 


15.44 ± 0.43 


1 


40 ± 0.26 


22.36 ± 1.66 


C Oph 


1 


,28 ± 0.25 


2 


.56 ± 0.14 


1 


.25 ± 0.14 


4 


23 ± 


13 





48 ± 


0.09 


6.28 ± 0.32 


0, 


36 ± 0, 


09 


9.28 ± 0.71 


0, 


17 ± 0.50 


< 18 


(T Ori 


1 


,56 ± 0.32 


2 


,71 ± 0,16 


0, 


,65 ± 0,12 


3 


50 ± 


17 





23 ± 


0.34 


< 6.7 


0, 


17 ± 


73 


< 13 








T ScO 


3, 


.28 ± 0.44 


3, 


,44 ± 0.15 


1 


,79 ± 0.14 


4, 


,73 ± 0, 


,13 


0, 


,84 ± 


0.09 


7.41 ± 0.26 


0, 


,58 it 0, 


,07 


10.62 ± 0.37 


0, 


,89 ± 0.38 


18.53 ± 2.98 


mean Thh^ 






2, 


.52 ± 0.44 






3, 


,60 ± 0, 


,68 








5.92 ± 0.94 








9.03 ± 1.19 






18.53 ± 2.98 


Note. — H/He rat: 


ios are ISM corrected and T/f 


H 


g are the derived X-ray temperatures in 


MK. The results for 


0^ Ori 


L C are not included 


in the mean T/f 


He ■ 



CO 



Table 6. Observed HEG H/He Line Ratios and Derived Thhs 



Star 


Neon 


Mag 


nesium 


Silicon 


Sulfur 




H/He 


'^H He. 


H/He 




H/He 


He. 


H/He 


^H He 


Tl range 




3.98 - 6.31 




^; Qi in nn 
O.ol - lU.UU 




10.00 - 15.85 




15.85 - 25.12 


Supevgiants (I, II) 


















C Pup 


1.21 ±0.16 


4.19 ±0.15 


0.26 ± U.Oo 


5.38 ± 0.22 


0.08 ± 0.04 


6.44 ± 0.59 


0.08 ± 0.29 


< 15 


Cyg \JdI JNo. 9 


















^^yg v^oz iNO. oa 






n Q/1 4- n oc; 
U.o4 ziz U.zo 


V QO -i- n T'i 

t .oJ ± v. 1 o 


0.94 ±0.19 


12.39 ± 0.90 


0.80 ±0.41 


17.66 ± 3.45 


Uri A 


0.68 ± 0.32 


3.50 ±0.50 


U.yO ± U.4d 


V iii J— 1 oo 


0.25 ±0.17 


8.13 ±1.62 


0.22 ±0.79 


< 20 




0.63 ± 0.10 


3.46 ±0.14 


U.zo zt U.Uo 


0.4y ± U.OO 


0.25 ± 0.16 


8.10 ±1.49 






e Ori 


0.53 it 0.11 


3.32 ±0.16 


A o^; _l_ nil 
0.26 ± U.ll 


0.32 ± U.56 


0.10 ±0.14 


< 8.4 






mean Thhh 




3.62 ±0.27 




6.24 it O.o4 




8.77 ± 1.61 




17.66 ± 3.45 




















HD150136 


0.67 ± 0.36 


3.46 ±0.56 


0.42 ± 0.10 


6.04 ± 0.41 


0.53 ±0.17 


10.27 ±0.99 


0.16 ±0.30 


< 16 


^ Per 


0.39 ± 0.12 


3.04 ±0.26 


0.17 ± 0.09 


4.81 ± 0.61 


0.26 ±0.27 


< 11 






t Ori 


0.80 ± 0.30 


3.68 ± 0.40 


0.64 ± 0.53 


6.39 ± 1.91 


0.39 ± 0.44 


< 12 






/3 Cru 


0.20 ±0.12 


2.53 ±0.38 


0.15 ±0.74 


< 7.6 










mean TuHe. 




3.18 ±0.46 




5.74 ± 0.57 




10.27 ±0.99 






Main Sequence 


















9 Sgr 


0.80 ± 0.34 


3.67 ±0.45 


0.34 ±0.27 


5.40 ± 1.33 


0.53 ± 0.68 


< 14 






HD206267 


0.93 ± 0.95 


< 5 


0.36 ±0.47 


< 7.5 


0.13 ±0.87 


< 13 






15 Mon 


0.40 ± 0.42 


< 3.8 


0.89 ± 2.23 


< 13 










6»i Ori 


2.92 ± 1.25 


5.68 ±1.06 


2.78 ±0.49 


11.68 ±0.93 


1.68 ±0.15 


15.43 ± 0.58 


1.00 ±0.19 


19.53 ± 1.44 


COph 


1.01 ±0.22 


3.95 ± 0.26 


0.67 ±0.17 


6.92 ± 0.51 


0.32 ±0.16 


8.82 ± 1.28 






a Ori 


0.37 ±0.13 


2.99 ± 0.28 


0.26 ± 0.39 


< 7 










T ScO 


1.88 ±0.31 


4.80 ± 0.28 


0.85 ±0.14 


7.44 ± 0.40 


0.67 ±0.11 


11.10 ±0.58 


0.15 ±0.24 


< 15 


mean Ttfife 




3.86 ± 0.60 




6.59 ± 1.01 




9.96 ± 1.21 







Note. — H/He ratios are ISM corrected and TnHe the derived X-ray temperatures in MK. The results for 6^ Ori C are not included 
in the mean Thhh- 



Table 7. Observed MEG and HEG G-Ratios 



Star 


Oxygen 


Neon 


Magnesium 


Silicon 


Sulfur 




MEG HEG 


MEG HEG 


MEG HEG 


MEG HEG 


MEG HEG 



Sit'pe.rgiants (I, II) 



Q Pup 


1 


10 


± 





18 





43 


± 





04 





65 ± 





07 





82 ± 





04 





78 


± 





06 





96 


± 





06 


1 


14 


± 





12 





83 


± 





49 


1.14 ± 0.42 


Cyg OB2 No. 9 

































49 ± 





18 












1 


95 


± 





46 
























Cyg OB2 No. 8a 

































40 ± 





07 





68 


± 





15 





89 


± 





08 


1 


32 


± 





21 


3 


16 


± 


1 


70 


1.98 ± 0.57 


S Ori A 





79 


± 





07 





60 


± 





07 





44 ± 





19 





62 ± 





06 





81 


± 





31 





90 


± 





14 





67 


± 





14 


2 


02 


± 


2 


19 


0.76 ± 0.66 


C Ori A 





96 


± 





06 





49 


± 





04 





67 ± 





06 





97 ± 





08 





78 


± 





10 


1 


22 


± 





16 


1 


08 


± 





26 


2 


60 


± 


5 


05 




e Ori 


1 


00 


± 





08 





86 


± 





06 





63 ± 





10 





75 ± 





07 





93 


± 





14 


1 


03 


± 





16 





43 


± 





09 


1 


60 


± 


1 


09 





Giants (IV, IB) 



HD160136 


1 


30 


± 





68 





87 


± 





16 





98 


± 





48 





65 


± 





10 





96 


± 





14 





76 


± 





08 





88 


± 





18 


1 


10 


± 





72 


0.6 


8 ± 


e Per 





52 


± 





12 





67 


± 





08 





73 


± 





18 





80 


± 





12 


1 


18 


± 





21 





92 


± 





14 


1 


52 


± 





99 





97 


± 





77 






L Ori 





91 


± 





10 





18 


± 





06 





32 


± 





09 





91 


± 





13 


3 


58 


± 


2 


20 


1 


86 


± 





61 


1 


05 


± 





37 
















(3 Cru 


1 


42 


± 





18 





92 


± 





12 





49 


± 





16 


1 


09 


± 





34 


3 


03 


± 


2 


56 


1 


31 


± 





66 



























Main Sequence 



9 Sgr 


1 


10 


± 





39 





44 


± 





08 





73 


± 





27 





25 


± 





09 





75 


± 





16 





54 


± 





12 





73 


± 





31 














HD206267 





20 


± 





09 





63 


± 





13 















66 


± 





12 


2 


57 


± 


2 


14 





67 


± 





19 


3 


66 


± 


3 


11 














15 Mon 





60 


± 





10 





88 


± 





17 





41 


± 





36 





53 


± 





16 


3 


33 


± 


4 


54 





60 


± 





26 





79 


± 





61 














el Ori C 


1 


34 


± 





78 





52 


± 





07 





89 


± 





37 





68 


± 





06 





92 


± 





15 





72 


± 





04 





92 


± 





07 





70 


± 





10 


0.84 ± 0.11 


C Opli 





93 


± 





16 





67 


± 





06 





51 


± 





09 





78 


± 





07 





98 


± 





16 





99 


± 





11 





S3 


± 





14 





36 


± 





21 




CT Ori 


1 


16 


± 





21 


1 


12 


± 





14 





56 


± 





14 





91 


± 





18 


1 


03 


± 





44 


1 


18 


± 





56 





74 


± 





36 














r Sco 


1 


57 


± 





20 





90 


± 





06 





90 


± 





13 





79 


± 





05 





57 


± 





06 





93 


± 





06 


1 


14 


± 





12 


1 


83 


± 





66 


0.64 ± 0.16 



a> 

Note. — G-Ratios are ISM corrected. ^ 
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Fig. 6. — MEG and HEG scatter plots showing the relation between Rfir and Rt=i for all 
luminosity classes. The error bars for Rt=i are determined for a X-ray continuum optical 
depth of 1 ± 0.5 and the adopted value of Rt=i is the average of these limits. Rfir data that 
only have lower limits or an uncertainty > lOi?* are not shown (see Tables 3 and 4). 
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Fig. 8. — MEG and HEG scatter plots showing the dependence of Tx (MK) (determined from 
H/He hne ratios) on the He-hke /ir-inferred radii, Rfir (determined from He-hke f /i hne 
ratios), for all OB luminosity classes. Rfir data that only have lower limits or an uncertainty 
> lOi?* and upper hmit Tx data are not shown (see Tables 3, 4, 5, and 6). 



